UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWh«n  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  report  number 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

Dynamic  Response  of  an  Ice-Breaker  Hull  to 

Ice  Breaking: 

5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

Final,  Aug.  1982  -  Feb.  1984 

«.  PERFORMING  ORG.  REPORT  NUMBER 

IIHR  Report  No.  273 

7.  AUTHORf*; 

A.  Mueller  and  R.  Ettema 

•.  CONTRACT  OR  GRANT  NUMBERf.J 

ONR:  683:ysh 

NR062-638 

N00014-79-C-041 1 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Institute  of  Hydraulic  Research 

The  University  of  Iowa 

Iowa  Citv.  Iowa  52242 

10.  PROGRAM  ELEMENT,  PROJECT,  TASK 

AREA  &  WORK  UNIT  NUMBERS 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Administrative  Contracting  Office 

12.  REPORT  OATE 

February  1984 

Office  of  Naval  Research 

536  South  Clark  Street,  Chicago,  ILL  60605 

13.  NUMBER  OF  PAGES 

114 

U.  MONITORING  AGENCY  NAME  A  ADDRESS^*/  different  from  Controlling  Office) 

15.  SECURITY  CLASS,  (of  thie  report) 

UNCLASSIFIED 

15a.  DECLASSIFICATION/ DOWNGRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thie  Report) 

Approved  for  Public  Release;  Distribution  Unlimited 

17.  DISTRIBUTION  STATEMENT  (of  the  abetted  entered  in  Block  20,  if  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

19.  KEY  WORDS  (Continue  on  revere*  eide  if  neceeeary  end  identify  by  block  number) 

Ice-Breaker  Ship,  Ship  Hull 
Interaction 

Motion,  Ice,  Ice  Breaking,  Hull-Ice 

20.  ABSTRACT  (Continue  on  revarme  aide  It  neceaaery  and  Identity  by  block  number) 

Experiments  were  conducted  with  the  objective  of  studying  the  dynamic  be¬ 
havior  of  an  ice-breaker  hull  moving  continuously  through  an  ice  sheet.  Of 
particular  interest  was  the  interaction  of  ice  forces  exerted  against  the  hull 
and  the  hull's  pitching  motion.  A  1: 48-scale  model  hull  of  the  USOGC  "Polar 
Star"  (WAGB10)  was  used  in  the  study.  The  model  hull  was  towed  by  a  towing 
carriage. 

cont . 

DD  ,  j°nM73  1473  ««T|ON_OM  NOV «  IS  OBSOLETE  UNCLASSIFIED 

S/N  0102-014*6601  .rnmiTV  CLASSIFICATION  OF  THIS  PAGE  (When  Dete  Entered) 


^CLASSIFIED _ 

i,k.CUWlTY  CLASSIFICATION  OF  THIS  PAGEfW7i»n  Data  Entmd) 

The  dynamic  response  of  the  model  hull  to  ice  breaking  was  used  to  infer 
the  vertical  ccmponent  of  the  ice  forces  exerted  against  the  hull;  to  evaluate 
the  temporal  variation  of  the  hull's  buoyancy  and  pitching  moment;  and,  to 
evaluate  the  inertia  force  components  of  the  hull's  motion.  Examination  of  the 
dynamic  response  of  the  model  hull  required  ccmputer-based  acquisition  of  tem¬ 
poral  records  of  its  pitch  and  draft  as  well  as  its  vertical  and  angular  ac¬ 
celerations.  Time  series  analyses  of  the  temporal  records  were  performed. 

It  was  found  that,  as  a  simplified  but  reasonable  approximation,  the  dy¬ 
namic  response  of  an  ice-breaker  hull  moving  at  constant  velocity  can  be  con¬ 
sidered  as  a  forced  oscillation.  It  is  characterized  by  the  ratio  of  the 
frequency  VQ/1C  to  the  hull's  natural  frequencies  of  pitching  and  heaving, 
which  for  the  scale  model  were  found  to  be  about  the  same  value,  fQ.  The 
frequency  VQ/1C  is  associated  with  the  cycle  of  ice  forces  that  is  experienced 
by  a  hull  moving  with  a  velocity  V0  through  a  sheet  of  ice  with  a  character¬ 
istic  length  lc.  The  mean  spacing  between  the  consecutive  cracks  developed  in 
the  ice  sheet  by  the  hull  is  proportional  to  lc;  consequently  the  parameter 
V0/lc  is  a  measure  of  the  dominant  frequency  of  ice  forces. 

For  relatively  low  velocities,  V0/f0lc«l,  the  force  terms  related  to  a 
hull's  buoyancy  are  largely  in  equilibrium  with  ice  forces  against  a  hull,  and 
the  inertia  forces  are  negligible.  If  V0/f0lc  =  1,  the  hull  is  in  a  state  of 
resonance.  As  the  buoyancy-force  components  are  180-degrees  out  of  phase  with 
the  inertia  force  components  and  the  two  sets  of  forces  partly  compensate  for 
each  other.  At  relatively  high  velocities  of  hull  motion,  Vo/f0lc>>l,  a 
superposition  of  two  conditions  occurs;  the  inertia  forces  compensate  the 
high-frequency  components  of  the  ice  forces  against  the  hull,  and  buoyancy 
forces  are  negligible;  and  transient  motions  of  the  hull  are  stimulated  at 
the  frequency  fG. 

The  bow  of  a  hull  moving  at  low  velocities  through  an  ice  sheet  depresses 
the  ice  sheet,  which  then  becomes  locally  flooded.  However,  for  a  hull  trav¬ 
elling  at  relatively  high  velocity,  there  may  not  be  sufficient  time  for  the 
ice  sheet  to  become  flooded.  Consequently,  the  hull  experienced  reduction  in 
buoyancy  in  the  region  of  its  bow,  and  an  increase  in  the  ice  forces.  The  re¬ 
duction  in  buoyancy  and  increase  of  vertical  ice-forces  were  estimated  from 
the  temporal  mean  values  of  forces  and  moments  acting  on  the  hull.  For  the 
model  hull,  it  was  estimated  that  the  depression  of  the  water  level  was  of 
the  order  of  0.01  m,  which  accounted  for  about  one  half  of  the  total  reduction 
in  the  hull's  displacement.  The  reduction  in  buoyancy  force  is  associated  witi 
both  increased  vertical  ice  forces  against  the  hull  and  increased  ice-hull 
friction. 
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i.  introduction 

The  use  of  physical  modelling  to  obtain  reliable  predictions  for  the  mean 
resistance  forces  that  are  encountered  by  an  ice-breaker  hull  moving  through 
ice  is  an  open  field  for  research.  An  improved  understanding  of  the  complex 
dependence  of  the  ice-related  resistance  force  Rj  on  the  parameters  of  the  ice 
sheet,  on  the  dimensions  of  the  hull  and  on  its  motion  requires  better 
knowledge  of  ice-breaking  mechanics.  This  report  presents  the  results  of  an 
experimental  study  in  which  more  parameters  were  measured  than  is  customary  in 
testing  practice.  One  of  the  study's  objectives  was  to  explore  whether  new 
insights  could  be  derived  from  such  additional  information.  In  particular, 
vertical  forces  due  to  ice  breaking,  buoyancy  and  inertia  as  well  as  the 
moments  involved  in  pitching  motion  of  the  hull  were  all  included  in  the 
analysis  of  its  dynamic  response  to  ice  breaking. 

Scaling  laws  provide  the  relations  with  which  to  scale  up  results  of 
model  tests  to  prototype  values.  Dimensional  analysis  is  the  standard 
approach  used  to  relate  the  resistance  force  Rj  to  key  parameters  involved  in 
the  motion  of  an  ice-breaker  hull  through  ice.  The  parameters  are  marshalled 
as  groups  of  dimensionless  variables  which  can  then  be  associated  with 
contributing  ice-resistance  forces.  This  approach  was  taken  by  Kashteljan 
(1968),  Lewis  and  Edwards  (1970),  Edwards  et  al.  (1972),  Enkvist  (1972),  Vance 
(1975)  and  others.  Edwards  et  al.  (1972),  for  example,  equated  the  resistance 
force  Rj  to  three,  ice-related,  resistance  forces,  which  they  determined  from 
the  results  of  model  and  full-scale  tests  of  the  USCGC  "Mackinaw"; 
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RI  =  CIB°£bh 


C__  p  bhV  / gh 
v  w  o  ^ 


+ 


P  gbh2 
w3 


resistance  due 
to  flexural 
failing  of 
ice  sheet 


resistance  due  to 
passage  of  the  hull 
through  broken  ice 
( velocity-dependent 
resistance) 


resistance  due  to 
interaction  with  ( 1 ) 

broken  ice  floes 
(gravity  resistance) 


...  2 

By  dividing  by  Pwgbh  ,  (1)  can  be  rewritten  in  the  following  nondimens ional 
form: 


*1 

pw9bh2 


=  C 


IB  p  gh 
wy 


+  C 


IS 


(2) 


where  p  =  density  of  water;  a  =  flexural  strength  of  ice;  b  =  beam  of  the 

W  £ 

hull;  g  =  acceleration  due  to  gravity;  h  =  ice  sheet  thickness;  VQ  =  velocity 
of  ship;  and,  C^g,  Cjg  and  Cy  are  coefficients. 

A  regression  analysis  can  be  used  to  provide  the  values  of  the  three 
coefficients  CIB,  Cjg  and  Cy  and  to  obtain  a  set  of  prototype-scale  data  from 
model-scale  tests.  The  three  coefficients  are  assumed  to  be  independent  of 
the  dimensionless  variables  in  (2).  For  a  more  detailed  discussion  of  the 
application  of  dimensional  analysis  to  the  study  of  the  resistance  to  motion 
of  an  ice-breaker  hull  moving  continuously  through  an  ice  sheet,  the  reader  is 
referred  to  the  aforementioned  studies  and  also  to  the  studies  by  Schwarz 
(1974),  Poznyak  (1981)  and  Enkvist  (1972,  1983).  The  latter  studies  also 
contain  descriptions  of  experiments  that  were  conducted  to  determine  the 
contributions  to  Rj  that  are  made  by  the  ice-related  resistance  components. 
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Resistance  to  hull  motion  through  level  ice  sheets  can  be  examined  more 
directly  by  analyzing  the  physics  involved  in  ice  breaking.  White  (1969) 
analyzed  vertical  and  horizontal  forces  and  the  moments  involved  in  the  ice- 
hull  interaction.  Enkvist  (1972)  gave  individual  estimates  of  the 
contributions  of  the  different  components  to  the  total  horizontal  resistance 
force,  Rj.  Milano  (1973,  1975)  tried  to  model  all  phases  of  the  ice-force 
cycle  for  a  hull  moving  through  a  level  ice  sheet  and  proposed  a  numerical 
model  for  predicting  RI#  The  present  study  follows  essentially  the  same  lines 
of  analysis  that  were  initiated  by  these  authors. 

An  ice-breaker  hull  moving  through  a  level  ice  sheet  breaks  ice  by 
applying  a  vertical  force  to  the  ice  sheet.  This  force  builds  up  to  the 
fracture  point  associated  with  the  flexural  strength  of  ice,  subsequently 
relaxes,  then  increases  until  the  ice  sheet  fails  once  again.  This  cyclic,  or 
time-dependent,  force  interacts  with  the  pitch  oscillation  of  the  hull.  In 
order  to  properly  understand  the  dynamic  behavior  of  an  ice-breaker  hull,  it 
is  necessary  to  consider  vertical  forces  and  the  pitch  moments  of  the  hull  as 
well  as  horizontal  forces.  In  the  present  study,  a  scale  model  of  an  ice¬ 
breaker  hull  was  instrumented  so  that  its  vertical,  horizontal  and  pitching 
motions  could  be  monitored.  This  was  accomplished  by  recording  the  horizontal 
force  imparted  to  the  hull;  the  draft  and  pitch  angle  of  the  hull;  and  the 
vertical  and  angular  accelerations  of  the  hull.  Fran  these  data,  the  behavior 
of  the  vertical  forces  imparted  to  ice  sheets  by  a  model  ice-breaker  hull  were 
estimated  as  a  function  of  hull  velocity.  Experiments  were  conducted  using  a 
1: 48-scale  model  hull  of  the  United  States  Coast  Guard  ice-breaker  ship  USOGC 
"Polar  Star"  (WAGB10). 
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II.  ANALYSIS  OF  FORCES  AND  MOMENTS  ON  AN  ICE-BREAKER  HULL 

The  following  simplified  analysis  of  forces  and  moments  on  an  ice-breaker 
hull  outlines  the  ideas  of  experimental  design  and  data  interpretation  that 
were  adopted  for  the  study.  The  analysis  is  applied  to  a  hull  moving  at 
constant  velocity  through  a  level  ice  sheet. 

Here  it  is  noted  that  differences  exist  between  the  dynamic  response  of 
hulls  which  are  towed  at  constant  velocity  and  self-propelled  hulls.  For  the 
latter,  velocity  may  not  be  constant  during  a  cycle  of  ice  forces.  Also,  the 
pitching  amplitude  experienced  by  a  self-propelled  hull  may  be  less  than  for  a 
hull  moving  with  constant  velocity  through  a  level  ice  sheet. 

Further,  because  the  ice-breaking  pattern  of  a  hull  may  be  influenced  by 
hull  shape,  so  the  ice- force  cycle  experienced  by  a  hull  may  be  affected  by 
hull  shape.  Rolling  oscillation  of  a  hull  may  also  affect  the  ice-breaking 
pattern,  or  cycle  of  ice  forces,  especially  when  the  natural  frequency  of  roll 
is  close  to  the  frequency  of  the  ice- force  cycle. 

An  additional  consideration  is  the  influence  of  the  material  behavior  of 
ice  on  the  cycle  of  ice  forces.  As  was  pointed  out  by  Milano  (1982),  the  cycle 
of  ice  forces  may  vary  for  changing  hull  velocity  because  the  deformation 
behavior  of  ice  may  change. 

A.  General  Considerations 

The  resistance  to  motion  of  an  ice-breaker  hull  is  the  horizontal 
component  of  the  total  force  R  due  to  hydranechanical,  viscous  and  ice  forces 
acting  on  the  hull.  The  total  force  is  the  integral  of  the  normal  and 
tangential  stresses  acting  on  the  surface  areas  of  the  hull,  and  can  be 
written  in  vector  notation  as 
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where  p^  =  hydrostatic  and  hydrodynamic  pressure;  p^  =  ice  pressure;  and,  dA  = 
elemental  area.  The  term  xv  dA  is  the  viscous  stress  which  is  negligibly 
small  compared  to  the  shear  stress  Xj  dA  attributable  to  friction  between  ice 
and  the  hull.  The  shear  stress  x  ^  can  be  expressed  as 


(2.2) 


where  V/|v|  is  the  unit  vector  in  the  direction  of  the  ice  movement  on  the 
surface  element  of  the  hull  and  F  is  the  coefficient  of  dynamic  friction 
which  may  be  dependent  on  the  relative  velocity,  V,  of  the  hull  and  ice. 

The  total  resistance  force  R  is  in  equilibrium  with  the  thrust  T  of  the 
propeller  as  well  as  with  the  weight  W  and  the  inertia  forces  of  the  hull 
(Figure  la).  The  equation  of  motion  of  the  hull  can  be  written,  in  vector 
notation,  as 


myXcg  =  R  +  W  +  T  (2.3) 

where  =  virtual  mass  of  the  moving  hull;  and  =  acceleration  of  the 
hull's  center  of  gravity.  The  corresponding  moments  of  the  forces 
R,  W,  and  T  have  also  to  be  in  equilibrium  with  the  net  moment  of  the  inertia 


force  acting  on  the  hull. 
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If  the  hull  is  floating  in  open  water  at  a  vertical  position  z  =  0  (see 
Figure  la),  the  hydrostatic  forces  B(z=0)  are  in  equilibrium  with  the  weight 
W  of  the  hull, 


W  +  B(z=0)  =  0  (2.4) 

If  the  bow  of  the  hull  rides  up  onto  the  ice  sheet  (Figure  16),  it  places  part 
of  the  hull's  weight  on  the  ice  sheet  and  the  hydrostatic  force  is  changed 
by  AB; 


B(z)  =  B ( z=0 )  +  AB  (2.5) 

The  ice  force  F  can  then  be  separated  from  the  total  resistance  force  R  and 
written  as 


F  =  R  -  B  =  //  p  dA  +  //  F  p  dA 

S  s  A  Ivl 


=  R  +  W  -  AB 


(2.6) 


so  that  the  vector  sum  of  total  force  R  and  weight  W  becomes 


R  +  W  =  Fjr  +  AB 


(2.7) 


Note  that  (2.6)  and  (2.7)  neglect  all  hydrodynamic  forces  associated  with 
viscous  drag  and  wave  resistance  to  hull  motion. 
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The  complex  spatial  and  temporal  distributions  of  pj,  fID  and  v  are  not 
known  and  cannot  be  measured  by  way  of  a  simple  experiment.  As  an  alternative 
approach,  considering  only  the  ice  forces  and  hull  motion  in  the  two- 
dimensional  x-z  plane  (see  Figure  lb),  the  ice  forces,  Fx,  Fz  and  their 
corresponding  moment  about  the  y-axis,  My,  can  be  estimated  from  the  dynamic 
response  of  the  hull  if  the  buoyancy  force  increment  AB  and  its  corresponding 
moment  increment  A ,  associated  with  pitching  and  heaving  of  the  hull,  can 
be  estimated  from  the  hull's  position  in  the  x-z  plane. 

B.  Estimation  of  Buoyancy  Forces 

The  equation  of  motion,  (2.3),  of  a  hull  moving  with  constant  horizontal 
velocity  along  the  x-axis  can  be  written  as 

"Vxcg  “  0  =  Fx  +  ABx  +  Tx 

and 

m  z  =  F  +  AB 
v  eg  z  z 

The  equation  for  the  pitching  motion  of  the  hull  is 

Iv6=My+ArVB  (2.10) 

where  Iv  =  virtual  mass  moment  of  inertia;  and  0  =  angular  acceleration.  The 

temporal  mean  values  of  z cg  and  0  are  zero,  since  heave  and  pitch  motions  are 
limited.  The  acceleration  x  is  also  zero  because  the  model  hull  is  towed  at 
constant  velocity. 


(2.8) 

(2.9) 
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The  hydrostatic  pressure  distribution  acting  against  the  hull  is  not 

disturbed  at  low  velocities  of  hull  motion  provided  the  depressed  portion  of 

the  ice  sheet  at  the  hull's  bcw  is  flooded  with  water.  At  higher  hull 

velocities  this  may  no  longer  be  true  because  there  may  not  be  sufficient  time 

for  water  to  flood  the  depressed  portion  of  the  ice  sheet.  In  addition,  it 

has  been  observed  in  field  and  model  tests  (e.g.,  Enkvist  1972)  that  a 

ventilated  area  may  form  between  a  high  velocity  bcw  and  the  depressed  portion 

of  an  ice  sheet.  The  disturbance  to  the  distribution  of  hydrostatic  pressure 

and  the  depression  of  the  water  surface  level  at  the  hull's  bcw  may  lead  to  a 

reduction  of  the  buoyancy  force  AB  and  give  rise  to  a  positive  value  for  the 

z 

buoyancy  force  component  AB^.  Consequently,  only  for  relatively  low  hull 
velocities  (see  Figure  1)  can  the  buoyancy  force  increments  AB  and  AB  ,  and 

Z  X 

their  resulting  moment  increment  AM^r,  be  estimated  using  the  following 
relationships : 


AB 


z 


A  p  gz 
w  w^3 


(2.11) 


and 


AB  =  0 
x 


(2.12) 


AMyB  =  "  Vw96 


(2.13) 


The  proportionality  constant  I  p  in  (2.13)  can  be  obtained  experimentally 
by  measuring  trim  angle  of  the  model  hull  due  to  a  known  weight  at  its  bcw  or 
stern.  Similarly,  one  can  determine  the  proportionality  constant  Aw  p^q  in 
(2.11). 
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where  Aw  =  area  of  the  waterplane  of  the  hull;  z  =  average  position  of  the 
hull  above  the  equilibrium  level  z  =  0;  Iyy  =  area  moment  of  inertia  of  the 
waterline  for  pitching;  and  0  =  pitch  angle.  The  additional'  terms  for  higher 
hull  velocities  can  only  be  crudely  estimated  since  no  measurements  are 
available.  The  additional  change  in  the  buoyancy  force  AB^  must  be  of  the 
order  of  the  reduction  of  the  weight  of  the  displaced  fluid;  that  is 

AB^  -  -  (2.14) 

where  Aj  and  A*z  are  the  characteristic  area  and  average  depth,  respectively, 
of  the  depressed  portion  of  the  ice  sheet.  The  moment  increment  AM^  can  be 
estimated  (Figure  lc)  as 

AM'  =  AB'L  (2.15) 

y  z 

where  L  =  the  moment  arm  of  A  B ' .  The  buoyancy  force  increment  A  B 1  is  of  the 

Z  X 

order  of  the  vertically-integrated  hydrostatic  pressure;  that  is 

t— 2 

AS'  ~  b  Pwg  (2.16) 

where  b  =  width  of  the  beam  at  the  relevant  station  of  the  hull.  The 
ratio  AB'/AB'  is  large  compared  with  unity  because  the  area  of  the  depressed 
portion  of  the  ice  sheet,  A-^ ,  is  large  compared  with  the  area  tA"z; 
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AB'  A,  A  z 
z  1 
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x  i  7 — 2 
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»  1 


bA  z 


(2.17) 


It  follows  from  the  preceeding  analysis,  (2.11)  through  (2.17),  that  the 
temporal  mean  values  of  Fx,  Fz  and  Mv  can  be  estimated  using  the  following 
equations : 


Fx  "  Tx  +  4Ex  <2.18) 

f2  =  +  V/5 

and 

M  =  I  p  g  0  +  AB'L 
y  yy  w  z 

The  force  ~F  is  approximately  equal  to  the  towing  force  T  because  the 
x  x 

contribution  of  the  buoyancy  term,  AB' ,  is  small.  The  force  F  is  increased 

X  z 

by  the  depressed  water  level  at  the  bow.  The  vertical  force  due  to  the  mean 

change  in  the  hull's  draft,  "z,  by  itself  is  an  underestimate  of  F  .  Part  of 

z 

the  buoyancy  which  had  compensated  for  the  hull's  weight  is  replaced  by  an  ice 
force.  Since  an  increase  in  ice  forces  against  the  hull  causes  a  larger 
friction  force  to  be  exerted  against  the  hull,  F^  and  T^  are,  therefore,  also 
increased  by  the  depression  of  the  water  level.  The  moment  M  is  under¬ 
estimated  when  it  is  equated  solely  to  the  term  containing  the  pitch 
angle  0 ,  if  the  water  level  of  the  bcw  is  depressed  at  higher  velocities.  For 
the  purpose  of  the  interpretation  of  the  experimental  results,  it  is  assumed 
that  the  distribution  of  hydrostatic  pressure  around  the  hull  is  undisturbed 
by  the  deflection  of  the  ice  sheet.  The  validity  of  this  assumption  is  to  be 


(2.19) 

(2.20) 
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checked  from  the  results  of  experimentation  and  estimates  of  AB'  AB'  and  AM' 

z  x  y 

are  attempted  by  making  use  of  relations  discussed  in  the  following  section. 


C.  Estimation  of  Ice  Forces 

The  components  Fx  and  Fz  of  the  resultant  ice  force  Fj  are  components  of 
the  same  distribution  of  the  ice  pressure  against  the  hull.  Because 
integration  of  (2.6)  for  each  force  component  involves  the  same  surface  area 
of  the  hull,  the  two  force  components  have  the  same  scale  and  their  ratio  can 
be  assumed  to  be  approximately  constant  for  each  hull  velocity, 

F 

=  o  (2.21) 

x 


This  assumption  is  valid  so  long  as  both  the  distribution  of  ice  pressure  and 
the  directions  of  ice-rubble  movement  around  the  hull  are  similar  from  test  to 
test. 

The  line  of  action  of  the  ice  force  usually  acts  through  the  hull's 
bow.  If  the  contribution  of  the  horizontal  ice  force,  Fx,  to  the  pitching 
moment.  My,  is  small  compared  with  the  contribution  of  the  vertical  ice  force 
Fz,  the  length  of  the  moment  arm  Lj  is  equal  to  the  distance  of  the  line  of 
action  of  Fz,  where  it  crosses  the  waterline,  to  the  hull's  center  of  gravity; 


L 


I 


(2.22) 


The  moment  arm  Lj  has  to  stay  within  one  half  of  the  ship  length. 
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The  vertical  ice  force,  Fz,  is  produced  by  the  buoyancy  of  the  depressed 
and  submerged  ice  and  the  bending  force  needed  to  break  the  ice  sheet; 

Fz  =  a2  Av/pw“pI*  9  h  +  a3a  h2  (2.23) 

where  and  are  constants  of  the  order  of  unity.  The  first  tern  in  (2.23) 
can  be  assumed  to  be  constant  with  time,  while  the  second  term  is  time- 
dependent  as  it  represents  the  increase  and  relaxation  of  the  bending  forces 
in  the  ice  sheet.  Inertial  forces  of  the  ice  floes  and  the  water  moving 
around  the  hull,  which  may  be  important  at  high  velocities  of  hull  motion,  are 
not  included  in  this  estimate. 


D.  Interaction  of  Ice  Forces  with  Heaving  and  Pitching  Motions  of  a  Hull 

The  equation  of  heaving  motion  for  a  hull  can  be  written  as 

it^z  +  c-jz  +  Pwg\z  =  Fz(t)  (2.24) 

(i)  (ii)  ( iii)  (iv) 

where  my  =  virtual  mass  for  vertical  motion;  =  area  of  the  waterplane  of 
the  ship;  and,  c^  =  damping  coefficient.  The  linearized  equation  for  pitching 
of  a  hull  can  be  written  as 


I  0  +  c J  +  p  gl  e  =  m  (t) 

v  2  w^1  yy  y' 

(i)  (ii)  (iii)  (iv) 

where  6  =  pitch  angle;  I  =  virtual  mass  moment  of  inertia;  I 

v  yy 

of  inertia  of  the  waterplane  about  the  y-axis;  c^  =  damping 


(2.25) 

area  moment 
coefficient; 
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and,  My( t)  =  moment  of  the  ice  forces  relative  to  center  of  gravity  of  the 
hull. 


Both  heaving  and  pitching  motions  of  a  hull  can  be  analyzed  as  forced 
oscillations.  The  inertia  terms,  (i),  the  damping  terms,  (ii),  and  the 
buoyancy  terms,  (iii),  are  in  equilibrium  with  forcing  functions,  (iv).  The 
natural  frequencies  of  heaving  and  pitching  of  a  hull  are  f^  and  fp, 
respectively,  where 


and 


1  /w 


0.5 


fh=7~<  * 
2n  v 


,  p  gi 

f  =  L-  (_w latj 

y  2ir  v 


0.5 


(2.26) 


(2.27) 


For  the  model  hull  of  the  USCGC  Polar  Star,  it  was  determined  that  f^  and  fp 
are  approximately  1Hz  and  are  both  hereinafter  referred  to  as  fQ.  The 
parameters  ir^,  Iv,  and  Iyy  are  properties  of  the  hull  and  were  estimated 
in  accordance  with  their  open-water  values.  The  damping  coefficients, 
c^  and  c^t  for  ice  are  as  yet  unknown.  However,  it  is  of  interest  to  examine 
the  form  of  a  solution  of  (2.24)  and  (2.25),  both  of  which  include  damping 
terms. 

With  a  forcing  function  in  the  form  of  Fz(t)  =  Fzo  exp(i2irft),  (2.24)  has 
the  following  steady-state  solution: 


z(t)  = 


zo 


4rr‘ 


,  2,,  2  ,2>2,  2,2. 

(mv  (fQ  -f  )  +  c,f  ) 


0^5  exp(i(2irft  -  a)) 


(2.28) 


where 
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2nf 

c 


(2.29) 


The  second  derivation  of  z(t)  with  respect  to  time  is 


z(t)  =  -  4ir2f2z(t) 


2 


(m  2(f  2-f2)2  +  c2f2) 

V  O  1 


q^5  exp(i(2irft  -  a)) 


(2.30) 


The  ratio  of  the  inertia  term  m  z(t)  and  the  buoyancy  term  p  gA  z(t)  is 

v  w  w 


m  z(t) 
v 


(2.31) 


p  gA  z(t) 
w3  w 


The  inertia  and  the  buoyancy  terms  are  180  degrees  out  of  phase  for  all 
frequencies  of  pitching  and  heaving  of  the  hull.  At  resonance,  f  equals  fQ, 
so  that  both  terms  have  the  same  amplitude.  At  low  frequencies,  when  f  tends 
to  zero,  the  buoyancy  term  pwgAwz  is  in  equilibrium  with  F(t)  and  the  inertia 


term  is  negligible.  At  high  frequencies,  when  f  tends  to  infinity,  mvz  tends 


to  F(t)  and  the  buoyancy  term  is  negligible.  The  phase  angle,  a,  of  z(t) 
relative  to  the  forcing  function,  F(t),  is  90  degrees  at  resonance.  It  tends 
to  zero  at  low  frequencies  and  to  180  degrees  at  high  frequencies,  for 
which  z(t)  is  in  phase  with  the  forcing  function. 

In  addition  to  the  above  steady-state  solution,  the  transient  solution 
which  is  the  response  to  an  impact  in  z,  or  z,  has  the  form 
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z ( t)  =  A  exp  ( 12  tt  (f^ - y — 2)  t)  (2.32) 

lv  16tt 

where  the  amplitude.  A,  depends  on  the  initial  conditions  of  z  or  z. 
Physically,  (2.32)  relates  to  a  damped  oscillation  at  a  fixed  frequency  which 
is  determined  by  the  resonant  frequency  fQ  together  with  the  ratio  of  the 
damping  coefficient  and  the  mass  mv*  The  ratio  of  the  inertia  and  the 
buoyancy  terns  is  close  to  unity  if  damping  is  small; 


mvz(t) 

Pw9V(t) 


,,2  2 
16ir  m 
v 


)  <  1 


(2.33) 


Analogous  relations  can  be  formulated  for  I  0  and  p  qi  0 . 

v  w  yy 

If  the  draft  z(t),  the  pitch  angle  0(t),  as  well  as  their  derivatives 

with  respect  to  time,  z(t)  and  0(t),  are  measured  as  functions  of  time, 

information  on  the  ice  forces  for  the  entire  of  frequency  band  can  be  gained 

by  Fourier  analysis.  A  direct  evaluation  of  the  sum  of  the  damping  term  and 

the  ice  forces  (terms  (ii)  and  (iv)  in  (2.24)  and  (2.25))  is,  however, 

••  •• 

difficult  because  the  phases  of  z(t)  and  z(t)  with  respect  to  6(t)  and  0  ( t ) 
must  be  measured  with  high  accuracy. 

The  temporal  behaviors  of  Iz(t)  and  My(t)  are  governed  by  the  breaking 
pattern  of  the  ice  sheet.  During  impact  with  the  ice  sheet,  the  hull  rides 
onto  the  ice  sheet  until  the  sheet  fails.  When  this  occurs,  the  vertical  ice 
force  acting  through  the  hull  relaxes.  The  radius  of  the  resulting 
approximately  circular  crack,  through  the  ice  sheet  in  the  vicinity  of  the 
hull's  bow,  scales  with  the  ice  sheet's  characteristic  length,  lc,  where 
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(2.34) 


and  Ef  =  elastic  modulus  for  flexure  of  the  ice  sheet;  v  -  EOisson's  ratio  for 
ice.  Hie  term  lc/v0  is  the  period  of  travel  between  consecutive  circular 
cracks.  Consequently,  the  time-dependent  part  of  the  forcing  functions  Fz(t) 
and  Mj,(t)  can  be  considered  as  a  series  of  impacts  of  different  intensities  at 
the  preferred  frequency  of  V!c.  Hie  dynamic  response  of  the  hull  to  ice 
breaking  is  thus  a  superposition  of  a  quasi-steady-state  solution  oscillating 

at  a  frequency  of  V(/lc  and  transient  solutions  starting  at  strong  impacts  and 

oscillating  at  a  frequency  of  ((f  ^  -  c^)/16tt 

o  1 

If  VVc  equals  uniby,  the  pitching  of  the  hull  is  in  resonance  with 
the  frequency  of  ice  breaking  a^  the  inertia  term  is  approximately  equal  to 
the  buoyancy  term  with  a  180-degree  phase  shift.  If  Vo/f0lc  is  significantly 
less  than  unity-i.e.,  for  lew  velocities-buoyancy  forces  are  in  equilibrium 
with  ice  forces  and  inertia  terms  can  be  neglected.  If  Vf0lc  is 

greater  than  unity-i.e.  for  high  velocities-inertia  forces  become  dominant 
and  are  chiefly  responsible  for  ice  breaking.  Draft  and  pitch  of  the  hull 
will  adjust  to  the  mean  forces  only,  provided  that  there  are  no  strong  impacts 
which  induce  a  superimposed  transient  solution  at  the  resonant  frequency. 

The  foregoing  analysis  of  the  dynamic  response  of  an  ice-breaker  hull  to 
ice  breaking  does  not  account  for  a  likely  feed-back  of  hull  motion  on  ice 
forces;  i.e. ,  in  the  foregoing  analysis,  the  z-axis  ccmponent  of  the  ice  force 

and  the  pitching  moment,  Fz(t)  and  My(t),  are  both  assumed  to  be  independent 
of  z(t)  and  6(f),  respectively. 


III.  EXPERIMENTATION 


A.  Experimental  Facilities 

1.  Hie  IIHR  Ice  Towing  Tank 

The  experiments  were  conducted  in  the  20m  x  5m  x  1.3-meter  deep  model  ice 
towing  tank  of  the  IIHR.  Figure  2  is  an  overview  of  the  towing  tank,  the  cold 
room  in  which  it  is  housed,  and  its  cooling  system.  The  cooling  system  is 
composed  of  two  compressors  which  provide  coolant  to  the  two  cooler  units 
situated  at  each  end  of  the  cold  roan.  The  compressors  are  in  turn  cooled  by 
water  pumped  fran  a  200  nr  sump.  If  the  sump  temperature  exceeds  a  certain 
limit,  a  cooling  tower  situated  outside  of  the  building  is  operated  to  cool 
the  sump  water. 

Fans  inside  the  four  cooler  units  draw  air  from  the  cold  room  and,  after 
the  heat  exchange  has  occurred  discharge  it  into  eight  ducts  which  extend  the 
whole  length  of  the  cold  roan.  The  chilled  air  is  forced  through  an  array  of 
half-inch  diameter  holes  along  the  base  of  each  duct,  thereby  producing  a  flow 
of  chilled  air  over  the  towing  tank.  The  four  ducts  of  each  side  alternate  to 
provide  an  even  distribution  of  the  cold  air.  Every  two  hours,  one  pair  of 
cooler  units  is  defrosted  by  electrical  heating.  Depending  on  the  ambient  air 
temperature  outside  the  cold  room,  the  total  cooling  capacity  of  the  system 
varies  between  15  and  20  kW,  and  enables  an  ice  sheet  to  grow  at  a  rate  of  1.5 
to  2  mm  per  hour. 

The  5-meter  wide  x  2.4-meter  long  motorized  carriage,  depicted  in  Figure 
3,  was  used  to  tow  the  model  ice-breaker  hull.  The  carriage  runs  along  rails 
on  the  concrete  walls  of  the  tewing  tank.  The  vertical  level  of  each  rail  was 
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adjusted  to  a  tolerance  of  ±  1.5mm  along  its  length.  An  angle  beam  on  one 
side  of  the  basin  gives  the  lateral  guidance  and  carries  the  rack  of  the  rack- 
and-pinion  drive  mechanism.  The  D.C.  motor  on  the  carriage  has  a  maximum 
torque  of  31  Nm  and  a  speed  range  of  58  to  1750  RPM.  A  1:15  gear  box 
increases  the  torque  to  413  Nm  and  gives  a  reduced  speed  range  of  3.9  to  117 
RPM.  The  effective  radius  of  the  pinion  is  0.06m;  consequently  the  carriage 
has  a  maximum  driving  force  of  6800  N  and  a  velocity  range  of  0.024  m/sec  to 
0.74  m/sec.  Higher  velocities,  up  to  2.2  m/sec,  can  be  achieved  if  a  1:5  gear 
box  is  coupled  to  the  D.C.  motor. 

In  order  to  measure  the  velocity  of  the  carriage,  a  wheel  carrying  a 
circular  array  of  holes  is  mounted  on  the  drive  shaft  of  the  D.C.  motor.  The 
passage  of  each  hole,  as  the  shaft  rotates,  is  sensed  by  a  photo  detector 
which  emits  a  light  through  the  hole.  The  number  of  pulses  counted  during  a 
time  interval  is  proportional  to  the  velocity  of  the  carriage.  The  length  of 
the  time  interval  is  0.371  seconds  so  that  1000  pulses  correspond  to  a 
velocity  of  0.333  m/sec.  After  each  interval  of  0.371  sec,  the  number  of 
pulses  is  latched  to  a  display  and  to  a  digital-analog  converter  which  holds 
the  voltage  during  the  following  interval  until  the  next  measurement  is 
available.  The  mean  velocity  of  the  preceding  interval  is  therefore, 
displayed  and  can  be  sampled. 

2.  The  Model  Ice-Breaker  Hull 

A  1: 48-scale  model  of  the  USOGC  "Polar  Star"  (WAGB10),  loaned  to  the  IIHR 
by  the  U.S.  Coast  Guard,  was  used  to  conduct  the  study.  The  abbreviated  lines 
of  the  prototype,  and  model,  hull  of  the  Polar  Star  are  depicted  in  Figure  4 
and  its  principal  model  dimensions  are  given  in  Table  1. 
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Care  was  taken  to  accurately  determine  the  relevant  mechanical  constants 
of  the  hull  in  order  to  accurately  evaluate  forces  and  moments  due  to  inertia 
and  buoyancy.  The  model  hull  was  loaded  with  a  mass  of  102.6  kg  so  that  its 
center  of  gravity  was  located,  to  within  a  tolerance  of  0.014  m,  at  the 
longitudinal  center  of  buoyancy.  The  area  A^  of  the  model  hull's  water  plane 
(the  plan  area  of  the  hull  at  its  water  line)  was  determined  to  be  0.81  m2. 
This  value  was  controlled  by  loading  the  floating  hull  with  a  mass 
increment  Am  and  measuring  the  change  in  draft  Az.  From  the  relation 


p  A  Az  =  Am 
w 


(3.1) 


the  result  =  0.82  m2  ±  0.08  m2,  was  obtained. 


The  hull  was  suspended  as  a  pendulum  so  that  its  moment  of  inertia  could 
be  evaluated  in  accordance  with  the  equation  of  motion 


(3.2) 


which  enables  I  to  be  determined  from  the  length  of  the  pendulum,  &  ,  the  mass 
of  the  hull,  m,  and  the  radial  frequency  of  oscillation,  to ,  where 


(3.3) 


The  length  H  Q  was  measured  by  loading  the  suspended  hull  with  a  mass  rrq  at  its 
bow  (see  Figure  5  for  a  definition  of  the  terms  in  (3.3)).  The  equilibrium  of 
moments , 
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m  A4  =  m^Jl  ^ 


(3.4) 


together  with  the  relation 

A4  =  04  (3.5) 

o 

permits  il  q  to  be  determined  in  terms  of  m,  rryi  and  the  angle  of 

rotation,  0.  The  resulting  moment  of  inertia,  I,  was  calculated  to  be  33.3  ± 

o 

0.5  kg  m  ,  corresponding  to  a  radius  of  gyration  k  =  0.57  m;  according  to  the 
relation 


I  =  mk2  (3.6) 

No  prototype  values  were  available  for  comparison  with  the  measured  values  of 
I  and  k. 

The  area  moment  of  inertia,  Iyy,  for  pitching  of  the  hull  about  the  y 
axis  (Figure  1)  was  determined  by  loading  the  floating  hull  with  a  mass  at  its 
stern  and  measuring  the  resulting  pitch  angle  of  the  hull.  A  value  of  Iyy  =■ 
0.25  ±  0.02  m^  was  calculated.  This  value  corresponds  to  a  metacentric  height 
GM  =  2.4  m.  The  vertical  distances  of  the  water-line  to  the  center  of  gravity 
and  the  center  of  buoyancy  were  assumed  to  be  small  compared  to  the 
metacentric  height  GM; 


GM  - 


(3.7) 
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The  added  mass  for  vertical  motions  of  the  hull  was  estimated  using  the 
method  formulated  by  Landweber  and  Macagno  (1957).  They  describe  the  hull  for 
each  section  as  a  Lewis  form  which  is  based  on  two  parameters;  the  ratio  of 
beam,  b,  to  the  draft,  H;  and  the  ratio  of  the  sectional  area  to  the  area 
defined  by  the  product  bH.  In  this  manner,  an  added-mass  coefficient,  a  = 
1.38,  was  calculated  for  the  model  hull.  The  added-mass  moment  of  inertia, 
Iv,  (or  the  virtual  moment  of  inertia)  was  estimated  by  measuring  the  pitch 
frequency  of  the  hull.  An  eccentric  driven  by  a  small  D.C.  motor  with  a 
variable  speed  control  was  mounted  at  the  stem  of  the  ship  hull  and  operated 
so  that  it  produced  an  oscillating  vertical  force  with  an  amplitude  of  about 
IN.  From  the  equation  of  motion  for  hull  oscillation,  (2.25),  the  resonant 
frequency  was  determined  to  occur  when 

1ir2f2  _  mg  GM  _  mg  GM  .  . 

p '  —  ram  (3,8) 

The  resonant  frequency,  fp,  was  determined  to  be  1.09  Hz,  which  is  associated 
with  an  added-mass  moment  of  inertia,  Iv  -  51.5  kg  m^,  or  an  added-mass  moment 
of  inertia  coefficient,  A  =  0.55.  The  two  coefficients,  a  and  A,  were 
estimated  for  hull  motion  in  open  water.  They  are  at  best  a  lower  limit  for 
the  coefficients  of  the  hull  when  it  is  breaking  an  ice  sheet. 

The  open-water  resistance  force  Fxo  was  measured  as  a  function  of  hull 
velocity.  The  resistance  force  Fxo  and  the  resistance  coefficient  C  are 
listed  in  Table  2  together  with  Froude  number,  Vy^gL,  and  Reynolds  number, 

V  L/v ,  for  each  test  run.  The  resistance  coefficient  C  was  estimated  as 
o 
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C  = - . (3.9) 

0.5  p  VS 
w  o 

where  S  is  the  wetted  area  of  the  hull.  As  was  assumed  for  the  analysis  of 
forces  on  an  ice-breaker  (Chapter  II),  it  is  evident  in  Table  2  that  the  open- 
water  resistance  to  the  hull's  motion  is  small  compared  with  the  ice-breaking 
resistance.  Furthermore,  the  Froude  numbers  of  the  hull  are  too  small  for  the 
wave  resistance  to  be  of  significance  (Lackenby  1965). 

The  coefficient  of  friction  for  evaluating  the  friction  force  acting 
between  ice  and  hull  was  measured  for  saline  ice  during  tests  conducted  by 
Arctec  Inc.  (Lecourt  and  Deslaurian  1976).  For  the  surface  finish  of  the 
model  hull  used  at  IIHR,  the  Arctec  Inc.  study  gave  FjD  =  0.478.  No 
additional  measurement  of  F^D  was  made  during  the  present  study. 

3.  Instrumentation 

The  model  hull  was  connected  by  a  25-millimeter  diameter  shaft  to  a 

dynamometer  supported  from  the  carriage,  as  is  shown  in  Figure  6.  The  shaft 

was  connected  to  a  vertical  plate  of  the  dynamometer  by  way  of  two  ball 

bushings  which  enabled  the  shaft,  and  the  hull,  to  move  vertically.  To  ensure 

an  almost  moment-free  connection  of  the  hull,  the  shaft  was  fixed  by  a  ball 

bearing  to  the  base  of  the  hull.  Friction  in  the  ball  bushings  affected  the 

vertical  motion  of  the  model  hull.  A  measurement  of  the  hysteresis  of  the 

vertical  position  of  the  hull  under  zero  load  resulted  in  an  uncertainty  Az  = 

±  0.65  mm  in  its  vertical  displacement.  This  value  corresponded  to  a  vertical 

resistance  force  A F  =  p  gA  Az  =  5N.  Once  static  friction  in  the  ball  bushing 

z  w  w 


23 


was  overcome,  the  vertical  resistance  is  reduced  because  friction  is  reduced 
to  a  lower  value  commens urate  with  dynamic  friction. 

The  dynamometer,  which  was  used  to  measure  the  translational  inertia 
force  of  the  hull,  or  the  towing  force,  consisted  of  two  0.15-millimeter  thick 
plates  which  allowed  a  horizontal  motion  that  was  restricted  by  a  50  lb 
Statham  UL4-50  load  cell  instrumented  with  a  Statham  UC-3  force  transducer. 
The  load  cell  limited  the  full-scale  displacement  of  the  dynamometer  to  0.12 
mm.  The  dynamometer  and  its  supports  had  a  resonant  frequency  of  6  Hz  which 
was  detected  in  the  signal  produced  by  the  force  transducer.  Yawing  of  the 
hull  was  restricted  by  the  location  of  a  vertical  rod  at  its  stem.  The  rod 
was  constrained  to  slide  in  a  36-millimeter  diameter  ring  which  was  connected 
by  a  frame  (which  can  be  seen  in  Figure  10 ) . 

Angular  and  vertical  accelerations  of  the  hull  were  measured  using 
accelerometers  mounted  at  its  bow  and  stem.  Each  accelerometer  was  formed  of 
a  5  lb  UL4-5  Statham  load  cell  loaded  with  a  mass  of  2.23  kg  in  order  to 
sensitize  it  to  vertical  acceleration.  The  range  of  the  force  transducers 
correspond  to  accelerations  of  about  lg,  or  9.81  m/sec^.  The  resolution  of 
each  accelerometer  was  of  the  order  of  0.01  m/sec  .  This  precision  could  not 
be  further  improved  since  measured  acceleration  amplitudes  were  small  compared 
to  the  maximum  load  range  of  lg.  Resonant  frequencies  of  the  accelerometers 
were  well  beyond  the  bandwidth  of  the  data  acquisition  system  that  was  used 
for  the  experiments.  The  accelerometers  were  sensitive  to  tilting  because  a 
roll  of  angle  <J>  would  reduce  the  vertical  load  against  the  hull  by  an  amount 

x 2 

A(gm)  =  gm(l-cos<|>)  =  gm 


(3.10) 
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This  reduction  of  the  load  is  equivalent  to  a  negative  (downward)  acceleration 
of 


9 


(3.11) 


The  accelerometer  resolution  of  0.01  m/sec^  is  equivalent  to  a  roll  angle 


of  <f>  =  2.6°.  In  seme  ice-breaking  tests,  a  lew  frequency  deviation  of  the 
vertical  acceleration  was  observed  and  was  attributed  to  rolling  of  the 
hull.  Both  accelerometers  sensed  the  same  roll  angle.  Therefore,  the 
measured  angular  accelerations  of  the  hull  were  not  affected  by  rolling  of  the 
hull.  Pitch  angles  were  too  small  to  be  sensed  by  the  accelerometers. 

The  force  transducer  of  the  dynamometer  and  the  two  accelerometers  were 
powered  by  a  bridge  amplifier,  which  was  located  outside  of  the  ice  room. 

Heave  and  pitch  motions  of  the  hull  were  measured  by  recording,  with 
linear  displacement  potentiometers,  the  vertical  position  of  the  hull  at  two 
positions.  The  potentiometers  were  excited  with  12  volts,  which  corresponded 
to  a  full  stroke  movement  range  of  0.15  m.  The  outputs  of  the  potentiometers 
were  transmitted  by  means  of  two  voltage  followers. 

The  six  voltages  from  the  load  cells,  accelerometers,  and  potentiometers, 
together  with  the  carriage  velocity,  were  scanned  using  a  digital  voltmeter. 
The  digitized  data  were  serially  transmitted  through  a  telephone  link  to  the 
IIHR's  HP-1000  computer  system  and  were  there  stored  on  disk.  The  bandwidth 
of  the  data  acquisition  link  was  100  Hz  so  that  each  channel  was  sampled  at  a 
rate  of  16.7  Hz.  At  this  sampling  rate,  the  6  Hz  resonance  on  the  force 
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signal  could  be  resolved  and  was  filtered  digitally  at  3  Hz.  The  signals  from 
the  two  accelerometers  were  filtered  with  a  48-dB/octave  analog  filter  with 
the  3-dB  point  set  at  6  Hz.  The  signals  of  both  the  linear  potentiometers  and 
the  carriage  velocity  did  not  need  filtering  before  sampling. 

B.  Experimental  Procedures 
1.  Ice  Sheet  Growth 

Ice  sheets  were  grown  from  a  1.3-percent,  by  weight,  urea  solution 
according  to  the  following  procedure:  With  the  cooler  system  operating  at 
full  capacity,  the  water  in  the  towing  tank  was  cooled  dcwn  to  a  temperature 
of  -0.2°C  and  the  roan  was  simultaneously  cooled  to  about  -12°C.  An  air- 
bubble  system  provided  the  necessary  mixing  of  the  solution  to  prevent 
supercooling  of  its  surface  layer.  Before  being  wet-seeded,  the  surface  of 
the  urea  solution  was  screened  to  remove  floating  ice  crystals.  Then,  the  air 
bubbles  and  the  blowers  of  the  cooling  units  were  shut  off  and  the  cold  room 
was  fogged  with  a  fine  spray  of  water  droplets.  The  spray  was  produced  with 
the  use  of  a  pressurized  air  spray  gun  and  a  pressurized  tank.  The  water 
droplets  froze  in  the  cold  air  and  settled  onto  the  surface  of  the  water  which 
had  by  then  reached  the  freezing  temperature  of  -0.4°C  for  the  urea 
solution.  The  wet-seeding  process  prevented  the  formation  of  relatively  large 
ice  crystals  and  enabled  a  multitude  of  small  crystals  to  grow  simultaneously 
over  the  surface  of  the  urea  solution. 

Each  ice  sheet  was  permitted  to  grew  to  about  80  percent  of  its 
thickness,  h.  The  roam  temperature  was  then  raised  to  about  2°C  and  the  ice 
sheet  allowed  to  warm  and  weaken.  The  bending  strength,  a _,  and  the  flexural 
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elastic  modulus,  Ef,  of  the  ice  sheet  were  monitored  until  af  reached  the 
values  prescribed  for  the  ice-breaker  test.  The  load,  f,  to  flexurally  fail  a 
cantilever  beam  of  length  l  and  width  d  was  used  as  to  estimate  the  value 
of  o^,  where 


°f  = 


(3.12) 


Four  to  six  cantilever  beams  were  failed  in  order  to  get  a  representative  mean 
value  of  Of.  The  flexural  elastic  modulus,  Ef,  was  determined  by  measuring  the 
increment,  6  ,  of  the  vertical  deflection  of  the  ice  sheet  due  to  small 
increments  of  a  point  load,  AP,  at  the  center-point  of  the  ice  sheet.  Thereby 


E 


f 


3_ 
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(1-v2) 

I  3 

pw9h 


(3.13) 


2.  Calibration  of  the  Transducers 

For  each  of  the  six  channels  of  data  (force,  accelerations  at  bow  and 
stem  of  the  hull,  two  linear  potentiometers  and  the  carriage  velocity)  the 
zero  level  and  the  sensitivity  of  the  transducer  were  determined  before  every 
run.  The  sensitivities  of  the  transducers  proved  to  be  constant  for  all  the 
runs. 

The  sensitivity  of  the  dynamometer  was  determined  to  be  33.5  N/Volt.  it 
was  measured  by  means  of  a  horizontal  load  applied  at  the  stem  of  the  hull. 
The  sensitivities  of  the  accelerometers  were  determined  by  measuring  the 
change  in  output  voltage,  AE,  due  to  an  applied  load  increment,  AF,  of 
1.47  N.  The  corresponding  sensitivity  of  -  0.42  m/sec^/Volt  equal  to  AF/m  AE 
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where  ma  is  the  suspended  mass  in  the  accelerometer.  The  sensitivity  of  the 
linear  potentiometers,  reckoned  to  be  13.0  mm/Volt,  was  measured  as  the  ratio 
of  a  change  in  length  to  a  change  in  signal  output.  The  sensitivity  of  the 
circuit  for  the  carriage  velocity  was  determined  by  correlating  its  output 
voltage  with  the  mean  velocity  of  the  carriage  (determined  by  use  of  a  stop 
watch  and  a  length  scale).  Then,  the  analog  setting  of  the  carriage  speed 
control  was  calibrated  against  the  output  voltage. 

Using  the  linear  displacement  potentiometers  to  measure  the  position  of 
the  model  hull  relative  to  the  carriage  rails,  it  was  found  that,  as  is  stated 
in  Section  III.A.l,  the  level  of  the  carriage  rails  varied  to  within  a 
tolerance  of  ±1.5  mm  along  the  length.  In  order  to  improve  the  accuracy  of 
the  measurements  of  hull  motion,  the  position  of  the  rails  relative  to  the 
water  surface  had  to  be  determined  as  a  function  of  the  position  of  the 
carriage  along  the  ice  tank.  This  was  accomplished  by  slowly  towing  the 
freely-floating  model  hull  in  open  water  and  observing  the  change  in  hull 
position  as  was  indicated  by  the  potentiometers. 

In  order  to  calculate  the  position  x(t)  of  the  carriage  during  the  ice¬ 
breaking  tests,  the  velocity  signal,  V(t),  was  integrated  with  respect  to 
time.  As  is  mentioned  in  Section  III. A. 3,  the  voltage  output,  Ev(t),  lagged 
behind  the  velocity  signal,  V(t),  because  Ev(t)  was  equal  to  the  mean  velocity 
of  the  preceeding  interval.  In  order  to  improve  the  estimate  of  the  mean 
actual  velocity  of  the  nth  time  interval,  the  step  in  the  signal  was  detected 
and  it  was  assumed  that  the  change  in  velocity  from  interval  n-1  to  n  was  the 
same  as  the  step  intervals  n-2  to  n-1.  An  accuracy  of  x(t)  of  ±  50  mm  was 
required  for  the  tests.  Therefore,  the  elevation  of  the  rails,  HR(x),  (Figure 
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7)  was  smoothed  by  a  moving-average  integration  technique  applied  to  50  mm 
lengths  of  rail,  so  as  to  be  consistent  with  the  resolution  of  hull 
position.  Table  3  is  a  summary  of  the  bandwiths,  and  the  sensitivities  of  the 
various  transducers  comprising  the  instrumentation.  The  errors  of  the  various 
measured  quantities  are  listed  in  Table  4. 

3.  Conduct  of  the  Experiments 

In  order  to  get  the  model  hull  prepared  for  each  test,  a  rectangular  slot 
was  cut  in  the  ice  sheet  while  it  was  warming  and  weakening.  With  the  hull  in 
open  water,  zero  levels  for  the  six  data  channels  were  recorded  and  the 
sensitivities  of  the  transducers  were  checked.  For  each  ice  sheet,  data  were 
taken  for  an  initial  run  of  about  two  hull  lengths,  with  the  hull  moving  at  a 
low  velocity  of  0.03  m/sec.  Subsequently,  a  second  run  was  executed  with  the 
velocity  set  as  is  indicated  in  the  list  of  parameters  given  in  Table  5.  In 
addition  to  recording  the  motion  of  the  hull  on  video  tape,  photographs  of  the 
hull  were  taken  during  each  test.  At  the  completion  of  a  test,  the  ice 
fragments  in  the  channel  behind  the  ship  were  also  photographed. 

4.  Numerical  Methods 

Mean,  H,  and  standard  deviation,  o  ,  of  a  time  series  x(t)  of  record 

X 

length  T  and  sample  rate  fd  were  calculated  using  the  relationships 


x  =  A  z.  x(n/£d) 

d  n=l 


(3.14) 


and 


2,0.5 


o 


(3.15) 
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The  covariance,  o 


xy 


of  two  time  series  x(t)  and  y(t)  was  determined  using 


°xy2  =  /t  £,  x(n/fd>  y(n/fa)  -  xy 


(3.16) 


d  n=l 


The  correlation  coefficient  for  two  time  series  x(t)  and  y(t)  is 


2 


o 


P 


xy 


(3.17) 


Power  spectra,  cross  spectra  and  coherence  functions  were  calculated  by  way  of 
the  following  methods  (e.g. ,  see  Bendat  and  Piersol  1971): 

The  time  series  were  split  into  N  segments  each  with  M  =  64  samples  or 
segment  period  Ts  =  M/f_j  =  3.84  seconds.  For  each  segment,  a  Fourier  series 
was  calculated  with  a  resolution  bandwith  Af  =  1/TS  =  0.26  sec.  Frequency 
leakage  was  reduced  by  subtracting  the  time  mean  and  by  a  cosine  taper  window. 

Each  tapered  segment  xs(t)  of  the  signal  was  then  represented  by 


n=-M/2  Ts 


M/2 

l 


Xs(nAf)  exp  (i2imAft) 


(3.18) 


where  X  (nAf)  is  the  Fourier  transform  at  nAf  and 
s 


T 

s 


Xg(nAf)  =  /  x  (t)  exp  (-i2irnAft)dt 


(3.19) 


o 


The  summation  over  the  range  [-M/2,  +  M/2]  is  sufficient  as  xg(t)  was 
filtered  at  f^/2  before  the  analog  to  digital  conversion  of  the  data.  The 


complex  number  (1/T  )  Xg  (nAf)  describes  the  amplitude  and  the  phase  of  the 
oscillator  at  the  frequency  nAf. 

The  one-sided  power  spectrum  Gx(f)  was  estimated  as 

2  N 

Gx(n&f)  =Nf~  Z  Xs(nAf)  xs(n&f)  (3.20) 

S  S=1 

where  n  =  1 ,  2 , . . . M/2 

and  G^(nAf)  is  the  density  of  the  power  in  the  frequency  space;  it  indicates 
the  contribution  of  the  frequency  interval  [(n-l/2)Af,  (n  +  1/2)  AF]  to  the 
mean  square  value 


- To  M/2 

(x(t)  -  x)  =  E  G  (nAf)  Af  (3.21) 

n=l  x 

The  one-sided  cross  spectrum  GXy  (nAf)  of  a  pair  of  functions  x(t)  and 
y(t)  was  estimated  by  calculating  the  Fourier  transform  of  segments  xs(t)  and 
ys(t)  and  using  the  relationship 

2  N  * 

G  (nAf)  =^|-  1  Xs  (rrtf)  Yg(nAf)  (3.22) 

y  S  S=1 

The  cross  spectrum  is  a  complex-valued  function  and  is  represented  as 


G  (nAf)  =  Co(nAf)  -  i  Quad  (nAf) 
xy 
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where  Co  (rAf)  is  the  cospectrum  and  Quad  (nAf)  is  the  quadspectrum  of  the 
functions.  The  magnitude,  G^CnAf)  is  equal  to  the  average  product  of  the 
magnitudes  lx  (nAf) I  and  |y  (nAf) I  and  the  phase  is  equal  to  the  average  phase 
difference  of  the  two  oscillators  at  (rAf).  The  covariance  is  equal  to  the  sum 
over  the  cospectrum,  or 


- - - -  m/2 

(x(t)  -  x)  (y(t)-y)  =  I  Co  (nAf)Af 

n=l  ** 


The  coherence  function 


is  defined  as 


(nSf) 


Gxy(tAf)2 

G  (rAf)  G  (rAf) 
x  y 


(3.23) 


(3.24) 


and  is  equal  to  unity,  if  the  phase  difference  of  the  two  oscillators 
Xs  (rAf)  and  Yg(rAf)  is  constant  for  all  segments. 
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IV.  RESULTS 


A.  Test  Parameters 

In  accordance  with  the  exploratory  nature  of  this  study,  only  a  limited 
set  of  data  was  taken.  The  hull  velocity  V0,  or  the  derived  frequency  scale 
V0/lc  (lc  =  characteristic  length  of  the  ice  sheet),  was  the  main  parameter 
that  was  varied  in  order  to  study  the  interaction  of  the  ice  forces 
experienced  by  the  hull  with  its  pitching  motion.  Two  runs,  or  tests,  were 
performed  for  each  ice  sheet:  The  first  run  involved  the  movement  of  the 
model  hull  at  the  relatively  low  velocity  of  0.03  m/sec;  and  a  second  run  was 
conducted  with  the  model  hull  moving  at  a  velocity  of  either  0.22,  0.3  or  0.6 
m/sec.  The  flexural  strengths  of  the  ice  sheets  were  chosen  so  that,  for  each 
hull  velocity,  runs  were  conducted  for  two  values  of  af.  Ice-sheet  thickness, 
h,  was  kept  constant  at  about  0.023  m,  with  the  exception  of  sheet  G  for  which 
h  =  0.032  m.  The  physical  parameters  of  the  seven  ice  sheets  used  in  the 
study  are  summarized  in  Table  5.  The  mean  and  standard  deviations  of  draft, 
vertical  acceleration,  pitch  and  angular  acceleration  of  the  model  hull  are 
presented  in  Appendix  1.  The  buoyancy  and  inertia  terms  are  additionally 
scaled  according  to  (2.24)  and  (2.25)  so  that  they  are  directly  comparable. 
The  scaling  factors  are  listed  in  Table  6.  Appendix  1  also  gives  the  mutual 
covariance  of  the  five  quantities  together  with  the  correlation 
coefficients.  The  errors  of  measurement  were  discussed  in  the  Section 
I I I. A. 3.  Table  4  is  a  summary  of  the  error  estimates. 

In  Appendix  2,  examples  of  the  time  series  of  the  measured  quantities  for 
every  run,  as  scaled  using  (2.24)  and  (2.25),  are  presented.  A  scaling  to 
prototype  values  was  not  the  purpose  of  this  study. 
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The  resultant  trends  that  were  obtained  for  the  temporal  mean  values  of 
forces  and  moments  experienced  by  the  hull  are  discussed  in  Section  B,  while 
the  time-dependent  characteristics  of  the  forces  and  moments  are  discussed  in 
Section  C. 

B.  Discussion  of  Mean  Forces  and  Moments 
1.  Calculated  Data 

The  measured*  temporal  mean  values  of  the  ice-related  resistance 

forces,  F  and  F  as  well  as  the  temporal  mean  values  of  the  moment  M  are 
xm  zm  c  ym 

represented,  in  Figures  8a,  b,  c,  as  functions  of  hull  velocity  VQ.  As  was 

discussed  earlier,  it  is  assumed  that  *F  is  approximately  equal  to  the  towing 

force  T^  (see  (2.18))  ;  that  "F  is  approximated  by  the  buoyancy 

term  p  gA  z  (see  (2.19));  and  that  M  is  approximated  by  the  buoyancy 
w  w  y 

term  I :  p^gS"  (equation  2.20). 

The  temporal  mean  horizontal  ice  force  F  generally  increases  with 

xm 

increasing  hull  speed.  The  measured  mean  buoyancy  force  F  increases  with 

zm 

hull  speed  up  to  VQ  =  0.3  m/sec;  and  subsequently  decreases  toward  VQ  =  0.6 

m/sec.  Fbr  an  ice-sheet  thickness  of  0.023  in,  the  moment  M  is  constant  for 

ym 

hull  speeds  up  to  a  value  of  0.3  m/sec  and  thereafter  increases  for  speeds 
upto  a  value,  VQ  =  0.6m/sec. 


* 


Measured  quantities  have  the  subscript  m. 
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2.  Comparison  with  Data  Obtained  by  Arc  tec  Inc. 

Extensive  model  tests  of  the  USCGC  "Polar  Star"  were  conducted  in  1976  by 
Arctec  Inc.  of  Maryland  USA  (Lecourt  and  Deslauriers,  1976).  For  these  tests, 
Arctec  used  saline  ice  (NaCl  solution),  instead  of  urea  ice,  with  flexural 
strength  ranging  from  7  to  12  kPa. 

Values  of  the  resistance  force,  F  as  were  obtained  in  the  present  study 
are  compared  in  Figure  9  with  those  obtained  during  the  tests  conducted  by 
Arctec.  The  IIHR  data  are  in  good  agreement  with  those  Arctec  data  that  are 
associated  with  a  coefficient  of  friction,  between  ice  and  hull,  of  F  = 
0.124.  The  Arctec  data  for  F  =  0.478,  which  was  associated  with  the  initial 
surface  finish  of  the  model  hull  (prior  to  its  use  by  the  IIHR)  are 
consistantly  higher  than  the  IIHR  data.  During  the  present  study,  no 
measurement  of  F  was  made  with  urea  ice. 

3.  Effects  of  the  Depression  of  the  Water  Level  in  the  Bow  Region 

The  trends  of  F  F  and  M  ,  that  are  shown  in  Figure  8,  can  be 
xni  zm  ym 

interpreted  by  considering  the  effect  that  the  depression  of  an  ice  sheet  has 
on  the  distribution  of  hydrostatic  pressure  that  is  exerted  against  the  hull 
(as  discussed  in  Section  II. B).  Because  little  information  is  available  on 
these  trends,  it  is  of  interest  to  get  at  least  approximate  estimates. 

At  the  relatively  low  value  of  hull  velocity  V0  =  0.03  m/sec,  the  loaded 
ice  sheet  was  flexed  down  below  the  water  level  and  became  flooded  with 
water.  The  photograph  Figure  10a  illustrates  the  flooding  of  such  a 
semicircular  area  which  had  a  diameter  of  about  twice  the  beam  of  the  hull. 
For  moderate  hull  velocities,  V0  =  0.22  m/sec,  this  area  did  not  flood  (Figure 
10b),  and  at  V0  =  0.6  m/sec,  not  only  did  flooding  not  occur  but  a  small  air 
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gap  was  visible  between  the  deflected  surface  of  ice  sheet  and  the  hull's  bow 

(Figure  10c).  It  can  be  concluded  from  these  observations  that,  at  V0  =  0.03 

m/sec,  the  distribution  of  the  hydrostatic  pressure  against  the  hull  is  not 

affected  by  the  depression  of  the  ice  sheet  and  that  the  analysis  as  is 

embodied  in  (2.7)  through  (2.12)  can  be  used  to  estimate  the  actual  ice  forces 

exerted  against  the  hull.  The  influence  of  water  surface  depression  due  to 

the  downward  flexing  of  the  ice  sheet  becomes  progressively  significant  at 

higher  velocities  of  hull  motion.  This  result  is  supported  by  the  observation 

that  at  Vn  =  0.03  m/sec,  the  measured  vertical  ice  force  F  correlates  well 

zm 

with  the  range  given  by  the  two  contributions  in  (2.23),  Fz  = 

2 

(pw-Pj)  gh  +  a f h  (Figure  11).  The  second  term  is  time  dependent  and 

contributes  only  with  its  temporal  mean  to  F~.  The  coefficients  a  and  a  in 

z  2  3 

(2.23)  were  set  equal  to  unity  for  this  comparison. 

In  order  to  analyze  the  effect  of  ice  sheet  depression  for  higher 
velocities  of  hull  motion,  estimates  of  the  ice  force  E'-  and  its  moment  M~  are 

z  y 

required  in  order  to  evaluate  AB'  and  AM'  (see  (2.18)  and  (2.19));  where 

z  y 


and 


F  =  F  +  AB' 
z  zm  z 


F 

zm 


+  V/ 


z 


M  =  M  +  AM'  =  M 

y  ym  y  y 


+  AB'L 
z 


The  estimates  can  be  derived  from  the  ratios  cm  =  f  /F  and  L  =  M  /F 

1  zm  xm  I  ym  zm 

wh i ch  are  plotted  in  Figures  12a,  b.  For  relatively  lew  hull  velocities, 
values  of  show  little  scatter  and  a  mean  value  of  a  ,  for  the  ice  sheets 
tested,  can  be  evaluated  as 
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al(Vo  =  0,03  ra/sec)  =  1-5  ±  °‘3  (4.1) 

It  is  evident  from  Figure  12  that  decreased  with  increasing  hull  velocity 
for  all  the  ice  sheets  that  were  tested.  However,  is  within  the  range  of 
values  given  in  (4.1),  except  when  V  =  0.6  m/sec,  for  which  the  value 
of  dropped  to  about  one  half  (sheet  F)  and  one  third  (sheet  E)  of  its  value 
given  in  (4.1). 

As  stated  in  Section  III.2.C,  it  can  be  assumed,  as  an  approximation, 
that  is  constant  provided  that  the  distribution  of  the  ice  pressure  and  the 
motion  of  the  broken  ice  is  similar.  If  it  is  further  assumed  that  =  1.5 
when  VQ  =  0.6  ra/sec  and  that  the  characteristic  area  A-^  of  the  depressed  water 
level  is  of  the  order  of  10-percent  of  Aw,  it  follows  from  (2.14)  and  (2.19) 
that 


AB'  =  1.5  F  -  F 
z  xm  zm 


(4.2) 


and 


1.5  F  -  F 
_ xm  zm 

°-l  V  wg 


(4.3) 


Estimated  values  of  AB^,  and  Az  are  listed  in  Table  7.  For  ice  sheets  E  and 
F,  the  estimated  value  of  Az  is  0.025  m.  This  value  seems  to  be  the  right 
order  of  magnitude.  It  is,  however,  a  quandary  as  to  how  the  ice  sheet  can 
transmit  the  total  vertical  force  Fz,  which  considerably  exceeds  its  bearing 
capacity  for  a  stationary  load. 
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Calculation  of  the  buoyancy  force  increment  AB^,  using  (2.16)  with  Az  = 
0.25m,  indicates  it  to  be  small  compared  to  the  measured  resistance,  or  the 
towing  force,  Tx.  The  approximation  F  =  T  is  therefore  a  valid  assumption 

A  X 

for  the  given  margin  of  experimental  errors. 


The  ratio  My/Fz  is  an  estimate  of  the  moment  arm  Lj  of  the  resultant  ice 

force  acting  against  the  hull.  When  VQ  =  0.03  m/sec,  the  mean  value  of  Lj  was 


Lj  =  0.73  ±  0.08  m 


(4.5) 


For  moderate  hull  velocities,  Lj  marginally  decreased  to  0.63  m  with  a 

root-mean  square  value  of  0.13  m.  When  V0  =  0.6  m,  Lj  increased  to  1.04m 

(sheet  F)  and  1.54  m  (sheet  E).  This  estimate  for  Lj  is  in  contradiction  to 

the  requirement  that  the  line  of  action  of  the  ice  forces  must  lie  within  the 

hull.  In  other  words,  this  result  supports  the  argument  that  F  must  under- 

zm 

estimate  the  true  value  of  the  vertical  ice  forces  acting  against  the  hull. 

The  moment  increment  AM^  can  be  estimated  using  an  argument  similar  to 
that  used  for  estimating  the  force  increment  A B' :  If  LT  =  0.73m  is  an  adequate 
estimate  of  the  moment  arm  for  all  hull  velocities,  then 


and 


AM'  =  I^a.,  F  -  M 
y  i  1  x  ym 


(4.6) 


The  moment  arm,  L,  of  the  buoyancy  force  increment  AB'  is 

z 


(4.7) 


38 


For  ice  sheets  E  and  Ff  L  =  0.35m  and  0.41  m,  respectively.  However,  these 
estimates  are  too  approximate  to  place  the  depression  of  the  ice  sheet  in  the 
bcw  region. 

The  results  summarized  in  Table  7  indicate  that  the  depression  of  the 
water  level  at  the  bcw  region  of  the  hull  depressing  an  ice  sheet  plays  an 
important  role  in  the  generation  of  ice  forces.  It  leads  to  an  increase  in 
the  resistance  to  motion  that  is  experienced  by  the  hull.  There  is  a  need  for 
more  detailed  experiments  to  investigate  the  influence  of  water  level 
depression,  at  the  bow  region,  on  the  dynamic  response  of  an  ice-breaker  hull 
to  ice  breaking. 


C.  Temporal  Variation  of  Forces  and  Moments 

In  section  II.D,  the  motion  for  the  heaving  and  pitching  motions  of  a 
hull  were  analyzed  as  forced  oscillations  and  the  following  relationships  were 
proposed : 


m  z  + 

V 

c^z  + 

p  gA  z 
w  w 

*  Fz(t) 

(2.24) 

10  + 

V 

c2e  + 

p  gl  6 
w3  yy 

-  y« 

Inertia 

Damping 

Buoyancy 

Forcing 

function 

(2.25) 

The 

vertical 

ice  force 

Fz(t)  and 

the  moment 

My(t) 

act 

as  forcing 

functions 

which  are 

in  equilibrium  with 

the  inertia 

forces 

and 

the  buoyancy 

forces. 


The  inertia  and  buoyancy  forces  were  measured  during  each  of  the 
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tests.  However ,  the  damping  terms  could  not  be  isolated  in  the  experiment. 
In  particular,  friction  in  the  ball  bushing  support  for  the  model  hull  caused 
a  fluctuating  vertical  force  of  unknown  size  to  act  through  the  center  of 
gravity  of  the  hull  and  to  damp  fluctuations  of  F(t).  Errors  due  to  this 
damping  effect  have  considerably  influenced  the  vertical  forces  that  were 
recorded  during  experimentation. 


1.  Amplitudes  of  Forces  and  Moments 

It  is  useful  to  compare  peak-to-peak  values  of  the  fluctuations  of  the 

inertia  and  the  buoyancy  forces  with  an  estimate  of  the  ice  forces  that  were 

needed  to  flexurally  fail  the  ice  sheet.  The  load  cycle  associated  with  the 

increase  and  subsequent  relaxation  of  the  ice  forces  was  considered  to  be  the 

source  of  the  temporal  fluctuations  of  the  forcing  functions  Fz(t)  and  M^(t) 

acting  on  the  hull.  In  Figures  13a  and  b,  typical  maxima  and  minima  of  peak- 

to-peak  values  of  PwgAwz(t)  and  myz  ( t ) ,  respectively  are  scaled  with  af h2,  and 

are  plotted  against  the  frequency  parameter  Vq/1c»  In  Figures  13c  and  d, 

typical  maxima  and  minima  of  peak-to-peak  values  of  pwgl  0(t)  and  1^6 (t) , 

respectively,  are  scaled  with  p  Ji2Lt  and  are  plotted  versus  V  A  .  The  moment 

L  o  c 

Lj  of  the  ice  force  acting  against  the  hull  is  taken  to  be  Lj  =  0.73  m. 

The  bearing  capacity  of  an  infinite  ice  sheet  sustaining  a  circular  load 
of  finite  diameter  is  reported  to  be  a  factor  of  1.5  to  2  times  afh  2.  For  a 
semiinfinite  ice  sheet,  the  factor  is  0.4  to  0.6  (Michel  1978).  The  buoyancy 
and  inertia  forces  that  were  measured  in  this  study  were  180  degrees  out  of 
phase  with  each  other.  Consequently  these  forces  partly  compensated  one 
another,  especially  at  the  resonant  frequency  (2.32).  The  observed  minimal 
peak-to-peak  values  appear  to  be  good  indicators  for  the  fracture  load  and  the 
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maximal  peak-to-peak  values  indicate  the  presence  of  a  strong  resonance.  From 

this  point  of  view,  the  fracture  load,  as  is  shown  in  Figure  13,  has  a 

2 

magnitude  between  0.4  and  0.8  times  a^h  ,  for  the  model  hull.  This  result  is 
in  agreement  with  the  fracture  load  of  a  semiinfinite  ice  sheet.  Strong 
resonance  in  the  pitching  motion  of  the  hull  was  recorded  for  values  of  Vq/1,^ 
between  1  and  2.  At  high  velocities  (V0/lc  =  3.4  and  4.2)  higher  peak-to-peak 
values  were  observed.  This  result  can  possibly  be  attributed  to  the  effects 
of  the  inertia  forces  of  the  ice  floes  and  the  water  underneath  the  ice. 

2.  Spectral  Distributions  of  Forces  and  Moments 

In  section  II. D  it  was  assumed  that  the  period  between  the  formation  of 
consecutive  cracks  in  the  ice  sheet  has  a  maximum  value  of  1C/V0 ,  where  lc  is 
the  characteristic  length  of  the  ice  sheet  and  VQ  is  the  speed  of  the  hull. 
If  the  rate  of  increase  of  the  ice  force  is  not  too  great,  the  spectrum  of  the 
forcing  function  has  a  peak  value  at  the  frequency  Vq/1^ 

The  response  of  the  hull  was  determined  by  calculating  the  power  spectra, 
the  cross  spectra  and  the  coherence  functions  for  the  temporal  records  of  the 
five  measured  variables  of  hull  motion — i.e. ,  draft,  pitch,  vertical  and 
angular  accelerations  and  towing  force.  The  power  spectra,  cross  spectra  and 
coherence  functions  for  the  measured  variables  are  summarized  in  Appendix  C. 

The  bandwidth  of  the  data  channels  are  listed  in  Table  3.  The  upper 
frequency  limit  of  Fx,  given  by  the  filter  frequency  of  3  Hz,  appeared  to  be 
inadequate,  especially  at  higher  velocities  (sheets  E,F)  and  for  sheets  of 
thicker  ice  (sheet  G).  The  draft  and  pitch  motions  of  the  hull  were  not 
filtered.  The  decline  of  their  spectra  is  therefore  given  by  the  motion  of 
the  hull.  The  accelerations  were  filtered  at  a  frequency  of  6  Hz.  The  low 
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frequencies  that  were  recorded  for  the  vertical  acceleration  z(t)  are  due  to 
rolling  motion  of  the  hull. 

Examples  of  the  time  histories  of  forces  and  moments  are  given  in 
Appendix  2.  Fran  these  examples,  the  relative  size  of  the  terms  in  (2.24)  and 
(2.25)  can  be  qualitatively  estimated.  The  power  spectra  (Figure  Cl)  give,  in 
a  more  quantitative  way,  the  average  contribution  of  each  frequency  interval 
to  the  mean  square  of  the  fluctuation  and  enable  a  comparison  to  be  made  of 

the  strength  of  the  forces  and  moments  in  the  frequency  space.  The  cross- 

•  •  •• 

spectra  and  coherence  functions  of  the  pairs  z(t)  and  z(t),  6(t)  and  9(t)  (see 
Figure  C2),  Fx(t)  and  z(t),  Fx(t)  and  z(t),  (see  Figure  C3),  Fx(t) 
and  9 ( t) ,  and  Fx  and  0(t)  (see  Figure  C4)  provide  information  on  the  frequency 
range  where  both  quantities  have  simultaneously  strong  spectral  lines,  and 
information  on  the  average  phase  relation.  The  cross  spectra  and  the 
coherence  functions  also  facilitate  the  checking  of  the  data  for 
consistency.  The  theoretical  coherence  functions  of  z  and  z,  as  well 
as  6  and  6 ,  are  unity.  Any  deviation  of  the  actual  values  from  the 
theoretical  values  must  be  due  either  to  unwanted  filtering,  to  noise  if  the 
signals  are  small,  to  error  of  measurements  or  to  a  violation  of  the  basic 
assumptions  which  lead  to  the  formulation  of  (2.23)  and  (2.24).  The  coherance 

function  of  0(t)  and  0(t)  (pitch  with  angular  acceleration)  is  close  to  unity 

in  a  frequency  band  around  the  resonant  frequency  fQ  =  1Hz.  The  coherance 
function  decreases  to  zero  for  higher  frequencies  of  pitching  motion.  The 
relationship  between  z(t)  and  z(t)  (draft  and  vertical  acceleration)  are  less 
coherent  if  the  fluctuations  of  the  forces  are  small.  This  may  be  due  to 

friction  in  the  ball  bushings  as  was  mentioned  before.  The  cross  spectra 

of  0(t)  and  0(t)  show  the  expected  behavior.  The  cospectrum  has  normally  a 
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negative  peak  at  the  resonance  frequency  and  the  quadspectrum  is  close  to 
zero.  This  result  agrees  with  the  fact  that  10  is  180  degrees  out  of  phase 
with  p  gl  6  . 

*  yy 

In  Appendix  C,  one  set  of  power  spectra  and  the  three  sets  of  cross 
spectra  and  coherence  function  are  organized  in  the  following  way:  Figure  (a) 
of  each  set  gives  examples  of  three  low  velocity  runs;  Figure  (b)  collects  the 
runs  at  medium  velocity  (VQ  =  0.22  m/sec  and  VQ  =  0.30  m/sec);  Figure  (c) 
shows  the  high  velocity  runs  (VQ  =  0.6  m/sec);  and  Figure  (d)  gives  the 
results  for  the  two  runs  of  sheet  G,  which  was  distinguished  by  its  greater 
ice  thickness. 

3.  Classification  of  Tests  According  to  Frequency  Response  of  the  Hull 

The  tests  can  be  classified  according  to  the  frequency  response  of  the 
hull  to  the  forcing  functions  exerted  against  it  during  the  process  of  ice 
breaking.  As  is  described  in  Section  II.D,  the  parameter  Vo/f0lc  can  be  used 
as  a  measure  of  the  dynamic,  or  frequency,  response  of  the  hull  to  ice 
forces.  It  was  found  for  the  model  hull  that  the  natural  frequencies  of 
heaving  and  pitching  are  sufficiently  equivalent  that  the  two  resonances 
cannot  be  separated.  Values  of  Vo/f0lc  for  the  tests  are  given  in  Table  5. 

For  the  relatively  low  hull  velocity  of  VQ  =  0.03  m/sec,  when  is 

significantly  less  than  unity,  the  acceleration  terms  are  small  and  the 
buoyancy  terms  are  in  equilibrium  with  the  forcing  function  due  to  ice 
forces.  The  power  spectra  of  the  inertia  terns  are  zero  and  the  power  spectra 
of  pitch  0  and  F^  are  nonzero  only  at  the  lowest  frequencies  (Figures 
Cla,d).  These  two  quantities  show  there  are  also  some  positive  values  in  the 
co-spectrum,  which  means  that  they  are  in  phase  (Figures  C4a  through  d). 
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An  average  distance.  Ax,  between  consecutive  cracks  in  the  ice  sheets  can 
be  estimated  from  the  fluctuations  in  the  temporal  records  of  pitch  angle  0 . 
Values  of  the  measured  frequency  of  the  ice  force  cycle,  f  =  Vq/Ax,  are 
compared  with  V0/lc  in  Figure  14.  In  Figure  15  the  same  data  is  used  to 
correlate  the  distance  x  directly  with  lc. 

At  hull  velocities  for  which  V0/f0lc  is  in  the  range  of  1  to  2,  buoyancy 
and  inertia  terms  are  of  equal  magnitude  and  partly  compensate  one  another,  as 
is  indicated  in  (2.33).  The  cross-spectra  and  power  spectra,  as  shown  in 
Figure  Cl-4b,  depict  the  resonance  effect,  especially  for  the  pitch  and 

angular  acceleration  of  the  hull  6(t)  and  0(t).  For  sheet  G  the  resonance  is 
also  evident  for  the  vertical  motion  of  the  hull,  z  and  z.  As  can  be  seen 
from  the  positive  resp.  negative  cospectrum  given  in  Figures  C4b  through  d 
(sheet  G) ,  Fx  is  in  phase  with  9(t)  and  z(t)  at  the  resonant  frequency  fQ  and 
out  of  phase  with  6(t). 

For  high  hull  velocities,  with  V0/f0lc  approximately  equal  to  4,  the 

power  spectra  of  the  buoyancy  terms  decrease  toward  higher  frequencies.  Ice 
sheets  are  broken  predominantly  by  inertia  forces,  which  show  strong  lines  in 
the  spectra  at  Vq/1c  (Figure  Clc).  For  sheet  F,  when  VQ  =  0.6  m/sec,  there  is 

a  peak  in  the  cospectrum  of  z  and  z  as  well  as  0  and  0  at  the  frequency  f  = 

1.3  Hz.  This  indicates  the  presence  of  transients  (see(2.32)  and  (2.33)). 
There  is  also  a  dominant  peak  in  the  power  spectra  of  Fx  at  the  frequency  f  = 
1.6  Hz  (Figure  Clc)  which  is  also  clearly  visible  in  the  time  histories 
(Figure  B).  This  peak  is  also  apparent  in  the  cross  spectra  of  Fx 
and  z ( t)  and  Fxand0(t)  (Figures  C3c  and  C4c).  The  phase  angle  is  not 
clearly  defined  by  the  ratio  of  the  cospectrum  and  the  quad-spectrum  because 
at  this  frequency  the  phase  error  of  the  data  acquisition  is  about  10  percent 
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of  2xt  .  The  distance  travelled  by  the  hull  during  the  period  T  =  0.62  sec, 
corresponding  to  the  dominant  frequency  of  ice  forces,  f  =  1.6  Hz,  is  equal 
0.36  m.  This  length,  which  is  2  times  lc,  is  an  indication  that  these 
fluctuation  of  Fx  may  be  produced  by  an  interaction  effect  of  the  pitch  motion 
and  the  breaking  of  the  ice  sheet.  As  was  stated  in  Section  II. D,  such  an 
interaction  effect  is  not  included  in  the  present  analytical  model  of  forced 
oscillation  of  an  ice-breaker  hull.  The  accuracy  of  the  spectra  is  marginal- 
limited  at  high  velocities  of  hull  motion,  because  the  observation  time  and 
data  acquisition  is  restricted  by  the  length  of  the  towing  tank. 

The  peaks  in  the  power  spectra  can  be  used  to  estimate  the  recorded 
frequency  f  of  the  forcing  function  and  to  estimate  the  distance,  Ax,  between 
consecutive  cracks  formed  in  the  ice  sheet  as  the  hull  moved  through  it.  For 
the  high  velocity  runs,  the  peaks  at  the  highest  frequency  were  chosen. 
Figure  14  shows  the  correlation  of  with  the  frequency  scale  Vo/lc  and 
Figure  15  shows  the  same  information  as  a  correlation  of  Vq/^  with  lc.  The 
good  correlation  that  is  evident  in  Figures  14  and  15  is  a  confirmation  that 
the  average  value  of  the  initial  breaking  length,  Ax,  is  adequately  estimated 
using  the  characteristic  length  of  the  ice  sheet,  lc.  This  result,  nonethe¬ 
less,  is  somewhat  contradicted  by  the  observation  that,  at  relatively  high 
hull  velocities,  the  mean  size  of  the  ice  floes  in  the  wake  of  the  hull  is 
significantly  smaller  than  lc.  Figure  16  shows  typical  sizes  of  these  ice 
rubble,  after  a  run  at  a  relatively  low  velocity,  VQ  =  0.03  m/sec  (Figure 
16a),  and  after  a  run  at  a  relatively  high  velocity,  V0  =  0.60  m/sec  (Figure 
16b).  However  subsequent  fracture  of  the  ice  floes  at  high  velocity  during 
passage  of  the  hull  may  explain  this  contradiction. 
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VI.  CONCLUSIONS 

Based  on  experimentation  involving  a  1: 48-scale  model  of  the  USOGC  ice¬ 
breaker  ship  "Polar  Star"  moving  with  constant  velocity  through  level  ice 
sheets,  the  following  principal  conclusions  were  drawn  on  the  dynamic  response 
of  an  ice-breaker  hull  to  ice  breaking: 

1.  The  dynamic  behavior  of  an  ice-breaker  hull  breaking  ice  can  be 
described  in  terms  of  a  forced  oscillation  and  classified  according 
to  the  parameter  V0/f0lc,  where  VQ  =  the  hull's  velocity,  fQ  =  the 
natural  frequency  of  pitching  (and  heaving)  motion  of  the  hull.  lc  = 
the  characteristic  length  of  the  ice  sheet,  relates  directly  to  the 
distance  between  consecutive  cracks  formed  circumferentially  around 
the  bow  of  an  ice-breaker  hull. 

The  flexural  failure  of  an  ice  sheet  by  a  hull  causes  a 
quasiperiodic  forcing  function  of  frequency  V0/ic  to  be  exerted 

against  the  hull.  The  amplitude  of  the  forcing  function  was 

.  .  ? 
estimated  to  be  in  the  range  of  0.4  to  0.8  a^h  . 

For  relatively  lew  hull  velocities,  VQ/f0lc  «  1,  the  buoyancy 
terns  in  the  equation  of  motion  are  approximately  equivalent  to  the 
ice  forces,  or  the  forcing  function,  exerted  on  the  hull  during  ice 
breaking.  When  an  ice  sheet  is  broken  by  a  hull  moving  at  relatively 
high  velocity,  V0/f0lc  »  1,  the  inertia  forces  of  the  hull  are 
principally  responsible  for  ice  breaking.  When  a  hull  is  moving  at  a 
velocity  such  that  V0/f0lc  =  1,  both  buoyancy  and  inertia  forces  of 
the  hull  are  large  compared  to  ice  forces.  These  forces  partly 
compensate  one  another,  because  they  are  180-degrees  out  of  phase. 
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2.  For  relatively  high  cruise  velocities,  the  hydrostatic  pressure 
distribution  acting  against  the  hull  may  be  reduced  because  the 
depressed  portion  of  the  ice  sheet  at  the  hull's  bow  may  no  longer  be 
flooded  with  water. 

The  distance  that  the  leading  edge  of  the  ice  sheet  was 
depressed  and  the  effect  on  the  hulls  motion  were  estimated  using  the 
assumption  that  the  ratio  of  vertical  and  horizontal  ice  forces  is  a 
constant  value  of  about  1.5.  The  depression  of  the  ice  sheet,  and 
the  water  level  was  estimated  to  be  of  the  order  of  0.01  m  for  a 
model  hull  moving  at  a  velocity  of  0.6  m/sec.  This  depression  of 
water  level  was  estimated  to  produce  about  one  half  of  the  reduction 
of  the  hulls  buoyancy  necessary  to  load  the  ice  sheet. 
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(a)  Total  forces 


(b)  Forces  and  moment  due  to  ice  breaking 


(c)  Influence  of  ice-sheet  and  water  level  depression,  at  the  bow, 
on  the  loading  of  a  hull. 

Figure  1.  Definition  sketch  of  forces  acting  on  an  ice-breaker 
hull. 
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Schematic  of  the  model-ice  towing  tank. 


(b)  Towing  tank  and  motorized  carriage 
Figure  2.  The  IIHR  model-ice  towing  tank 
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Figure  4.  Abbreviated  lines  of  the  USGSC  ice- 
breaker  ship  ’Tolar  Star” 
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Figure  5.  Definition  sketch  of  the  pendulum  arrangement 
that  was  used  to  determine  the  mass  mcment  of 
inertia  of  the  model  hull. 
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Figure  6.  Details  of  the  dynamometer  that  was  used  for  measuring 
towing  force. 
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Figure  7.  Relative  height,  Hp(x) ,  or  vertical  correction  for  position  along 
the  tracks  of  the  motorized  carriage. 
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Figure  8.  Temporal  mean  values  of  measured  forces  and  moments  acting  against  the  hull. 
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Figure  9.  Comparison  of  F^  with  the  results  of  the  model  tests  conducted  by  Arctec  Inc.  The  number 
on  each  plot  point  indicates  ice  sheet  thickness  in  mm. 
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(a)  For  a  low  speed 
(V  =  0.03  m/s),  the 
depression  is  flooded 


(b)  For  a  medium  speed 
(V  =  0.22  m/s),  the 
depressed  area  is  dry 


(c)  For  a  high  speed 
(V  =  0.60  m/s),  an 
air  gap  was  visible 
at  the  bow 


Figure  10.  Depression  of  the 
ice  sheet  at  the  bow  of  the 
hull. 


.  Comparison  of  P^,  at  VQ  =  0.03  m/s, 
with  Aw(pw  -  P;[)  gh  +  of h2. 


Figure  11 


(a)  Force  ratio  versus  hull  velocity 


Hull  Velocity,  V0  ~ (m/sec) 


(b)  Moment  aim  versus  hull  velocity 


Figure  12 .  Variation  of  force  ratio  «f  and  moment  arm  LT  with 
hull  velocity,  V  .  1  1 
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Figure  13.  Typical  maxima  (0)  and  minima  (A)  of  peak-to-peak 

values  of  vertical  forces  and  moments  plotted  versus  frequency 
scale  of  ice  breaking,  V  /I  . 
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Frequency  of  Ice  Breaking,  V0//c~(Hz) 

Figure  14.  Correlation  of  frequency  of  forcing  function,  f  , 
with  frequency-scale  of  the  breaking,  Vq/1  .  m 


Characteristic  Length,  2C 


Correlation  of  ice-breaking  length.  Ax,  with  character¬ 
istic  length  of  ice  sheet,  1  . 

c 


Figure  15. 
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(b)  Hull  moving  at  a  relatively  high 
velocity  Vq  =  0.60  m/s 


Figure  16.  Typical  sizes  of  ice  rubble 
after  the  model  hull  has 
transmitted  ice  sheets. 
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Table  1.  List  of  Physical  Parameters  of  Ship  Model 


Overall  length  of  hull 

2.534  m 

Length  between  perpendicular 

2.235  m 

Beam  (maximum  width) 

0.530  m 

Beam  (design  at  waterline) 

0.495  m 

Draft  (design) 

0.178  m 

Displacement 

m  =  102.6  kg 

Virtual  mass  coefficient  (heaving  motion) 

a  =  1.38 

Mass  moment  of  inertia 

I  =  33.3  kgrtV 

Virtual  mass  moment  of  inertia  coefficient 

A  =  0.55 

Area  of  the  waterplane  of  the  hull 

Aw  =  0.82  m2 

Area  moment  of  inertia  for  pitch  motion 

Iyy  =  0-25  ra 

Wetted  area  of  hull 

S  =  1.32  m2 

Table  2.  Measurement  of  Open-1 


Hull 

Velocity 

Vo 

(m/sec) 


0.03 

0.0064 

37 

X 

103 

0.15 

0.032 

188 

X 

103 

0.22 

0.046 

276 

X 

103 

0.30 

0.064 

377 

X 

10  3 

0.59 

0.128 

741 

X 

103 
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Resistance  to  Hull  Motion 


Resistance  Coefficient  of 
Force  Resistance 


Fx 

(N) 

C 

0.2 

0.34 

0.4 

0.023 

0.4 

0.012 

0.5 

0.0084 

1.5 

0.0065 
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Table  3.  List  of  the  Bandwidths  and  the  Sensitivities 


Bandwidth 

Range  Limited  by 

(Hz) 


Force  Fx  0-3 

Accelerometers  0-6 

z(t)  0.5-6 

Linear  0  -  00 


Potentiometers 


digital  filtering  to 
remove  resonant 
frequency  and  6  Hz 

analog-filter  before 
A/D  conversion 

Rolling  of  hull 


Speed  control  of  the 
carriage 


of  the  Instruments 


Sensitivity 


33.5  N/Volt 

0.42  m/sec ^/volt 


0.013  m/Volt 


Velocity 


0 


0.137  m/sec/Volt 
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Table  4.  List  of  Experimental  Errors 


Parameter 


Error 


Main  Source  of  Error 


rx 

Draft  z 

pw9V 

Pitch  0 
p  gl  0 

w  yy 


mvz 


±  IN 

<  ±0.00065  m 
±  5N 


electrical  drift 

static  friction  of 
ball  bushings 


±  0.0004  radius 
+  iNm 


electrical  drift 


±  0.02  n^/sec^  Rolling  of  Ship 

••  o 

±  5N  (Az  =  0.005  m/sec  /degree 

error  occuring  at  low  frequencies) 

±  0.01  sec  ^  electrical  drift 


±  0,5  Nm 
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Table  5.  Physical  Parameters  of  Ice  Sheets  and  Test  Runs 


Sheet 

Date 

h 

°f 

(kPa) 

Ef 

lr 

ofh2 

Run  No 

v_ 

V 

o 

V 

o 

(m) 

JL 

(MPa) 

c 

(m) 

f 

(N) 

lc 

(m/sec)(sec  -*-) 

f  1 
o  c 

A 

10/15/82 

0.25 

23 

11 

0.20 

14.4 

02015 

03015 

0.03 

0.22 

0.15 

1.1 

0.15 

1.1 

B 

10/27/82 

0.021 

29 

6.5 

0.15 

12.9 

04027 

06027 

0.03 

0.22 

0.20 

1.5 

0.20 

1.4 

C 

10/18/82 

0.023 

10 

4.2 

0.15 

5.3 

02018 

03018 

0.03 

0.30 

0.20 

2.0 

0.20 

2.0 

D 

10/25/82 

0.023 

25 

6.0 

0.16 

13.0 

04025 

05025 

0.03 

0.30 

0.19 

1.9 

0.19 

1.9 

E 

10/20/82 

0.022 

19 

3.8 

0.14 

9.3 

04020 

05020 

0.03 

0.59 

0.21 

4.2 

0.21 

4.2 

F 

10/29/82 

0.023 

29 

8.5 

0.18 

15.3 

06029 

07029 

0.03 

0.60 

0.17 

3.4 

0.17 

3.4 

G 

10/22/82 

0.032 

18 

8.3 

0.22 

18.1 

04022 

05022 

0.03 

0.22 

0.13 

1.0 

0.13 

1.0 

T 

(sec) 


47.7 
36.0 

48.9 
19.2 

42.9 

29.8 

44.1 
30.0 

59.1 

6.0 

49.1 
12.0 
49.1 
27.0 


Table  6.  List  of  Scaling  Factors 


Draft: 

Vertical  Acceleration: 
Pitch: 


P_.g  A  z  =  7950  z  (N) 

w  w 
••  •« 

n^z  =  244  z  (N) 

P  gl  0  =  2450  9  (Nm) 

w  yy 

Iv0  =  51.6  0  (Nm) 


Angular  Acceleration: 
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Table  7.  Effect  of  the  Depression  of  the  Water  Level  at  the  Bow 


Sheet 

Vo 

Fx 

F 

z 

M 

y 

al 

LI 

AB' 

z 

Az 

AM' 

y 

L 

AB' 

X 

(m/sec) 

(N) 

(N) 

(Nm) 

(m) 

(N) 

(m) 

(Nm) 

(m) 

(N) 

E 

0.59 

18.6 

8.74 

13.72 

0.47 

1.56 

19.2 

0.024 

6.7 

0.35 

1.5 

F 

0.60 

26.5 

19.9 

20.8 

0.75 

1.04 

19.9 

0.025 

8.2 

0.41 

1.5 

APPENDIX  A 


Mean,  Standard  Deviations,  Covariances  and  Correlation  Coefficients 

of  the  Experimental  Runs 


Nomenclature : 
Draft 

Az 

Pitch 


z(t) 

z(t) 

0(t) 

e  (t) 


Omega 


71 


SHEET  A 

VEL  = 

0.03 

(RUN  NO 

:  C2015 

) 

F  x 

Draft 

Az 

Pitch 

Omega 

N 

m 

m/sec**? 

radians 

1/sec  *  *  2 

KEAN 

13.64 

0.0076 

-0.007 

0.0043 

-n.007 

SIGMA 

3.41 

0.0 0 E  2 

0.003 

0.0009 

0.C04 

scaled  as 

forces  and 

moments 

N 

N 

N 

Nm 

Nm 

MEAN 

1 

7.6  4 

20. r9 

-n.8i 

10.46 

-0.38 

SIGMA 

3.41 

1.27 

0.71 

2.31 

0.21 

Covariance 

of 

F  *  with 

:  <N*hi)  or 

( N*  Nm ) 

-0.58 

-C.83 

4.25 

-0.09 

Draft  with 

:  (N+N )  or 

(  N  +  N  m  ) 

-0.20 

1.60 

-0.07 

A  z  with 

:  CN*Nm) 

-0.79 

0.02 

Pitch  with 

;  ( N  m  *  n  ro ) 

-0.17 

Corr.  Coefficient 

F  x  with 

t 

• 

-0.13 

-0.35 

0.55 

-0.17 

Draft  with 

• 

• 

-C  .23 

0.54 

-0.26 

A  z  with 

• 

• 

-0.46 

0.15 

Pitch  with 

♦ 

-0.35 

SHEET  A 

VEL  = 

0.22 

(RUN  NO 

:  03015 

) 

Fx 

Draft 

Az 

Pitch 

Omega 

N 

m 

m / sec**2 

radians 

1/sec**2 

MEAN 

13.59 

0.0025 

-0.01? 

0.0050 

-0.008 

SIGMA 

3.76 

0.0004 

0.01  0 

0.0017 

0.050 

scaled  as 

forces  and 

moments 

N 

N 

N 

Nm 

Nm 

MEAN 

13.58 

20.19 

-2.92 

12.27 

-0.40 

SIGMA 

3.76 

3.15 

2.43 

4.23 

2.57 

Covariance 

of 

F  x  with 

:  (N*N)  or  ( N*Nm ) 

7.16 

-0.86 

9.19 

-2.77 

Draft  with 

1  (h*N)  or  (N*Nm) 

0.08 

6.48 

-1.71 

A  z  with 

:  (N+Nm) 

0.43 

0.46 

Pitch  with 

:  ( N  m  *  N  m ) 

-7.87 

Corr.  Coefficient 

F  x  with 

• 

« 

0.60 

-C.08 

0.58 

-0.29 

Draft  with 

* 

0.01 

0.64 

-0.21 

A  z  with 

* 

♦ 

0.G4 

0.07 

Pitch  with 

• 

-0.72 

SHEET  B 

VEL  =  0.03 

(RUN  NO 

I  04  Q?7 

) 

Fx 

Draft 

A  z 

Pitch 

Omega 

K: 

m 

m/sec**? 

radi ans 

1  / s  ec  *  *  2 

KEAN 

3.06 

n.aoi7 

0 .00A 

0.0046 

-0.004 

SIGMA 

1.65 

0.000? 

0 .005 

0 .0008 

C.01 1 

scaled  as 

forces  end  moments 

N 

N 

N 

Nm 

Nm 

KEAN 

8.06 

13.12 

0.92 

11.20 

-0.23 

SIGMA 

1  .65 

1  .7? 

1.24 

2.03 

0.5  8 

C  ovariance 

of 

F  x  with 

I  (  N  *  N )  or  ( N  *  N  m ) 

0.52 

-0.47 

1  .  32 

-0.05 

Draft  with 

;  ( N  *N  )  or  ( N  *  Nm ) 

-1.14 

1 .  PC 

-0.16 

A  z  with 

:  (N  *  N  rr ) 

-1  .30 

0.06 

Pitch  with 

:  (Nm*Nm) 

-0.53 

Corr.  Coefficient 

F  x  with 

♦ 

• 

0.24 

-0.23 

0.39 

-0.05 

Draft  with 

• 

• 

-0.70 

0.67 

-0.21 

A  z  with 

-0.52 

0.08 

Pitch  with 

l 

-0.4  5 

SHEET  B 

VEL  = 

C  .22 

(RUN  NO 

:  06027 

) 

Fx 

Draft 

A  z 

Pi  t  ch 

Omega 

M 
•  * 

m 

m/sec**2 

radians 

1 /sec**? 

MEAN 

15.1? 

C .  0  0  2  7 

-0.001 

0.0047 

-0.001 

SIGMA 

3.60 

0.0003 

0.008 

0.0014 

0.059 

scaled  as 

forces  and 

moments 

N 

N 

N 

Nm 

Nm 

MEAN 

1 

5.1? 

21.51 

-0.2  6 

11.64 

-0.04 

SIGMA 

3.60 

2.61 

1  .91 

3.33 

3.05 

Covariance 

cf 

F  x  with 

:  <N*N>  or 

(N*Nm) 

0.62 

-2.43 

5.17 

-2.77 

Draft  with 

:  (N+N)  or 

( N  +  Nm) 

-0.91 

2.67 

-1  .66 

A  z  with 

:  (  N  *  H  m  ) 

-2.73 

3.00 

Pitch  with 

:  (  n  m  *  n  m ) 

-7.39 

Corr.  Coefficient 

F  x  with 

t 

♦ 

0.07 

-0.35 

0.  A3 

-0.2  5 

Draft  with 

# 

* 

-  u .  1  a 

0.31 

-0.21 

A  z  with 

• 

• 

-0.43 

0.51 

Pitch  with 

» 

• 

-0.73 

SHEET  C 

VEL  = 

G  .03 

(RUN  NO 

:  02018 

) 

Fx 

Draft 

Az 

Pitch 

Omega 

N 

m 

m/se  c**2 

radians 

1 /sec**2 

MEAN 

7.73 

o.oou 

0.004 

0.0038 

-0.002 

SIGMA 

2.01 

0.0001 

0.003 

0.0006 

0.G04 

scaled  as 

forces  and 

moments 

N 

N 

N 

Nm 

N  m 

MEAN 

7.73 

10. 94 

1.06 

9.21 

-0.12 

SIGMA 

2.01 

0.66 

0.68 

1.35 

0.23 

C  ovar iance 

of 

F  x  with 

:  (N*N )  or 

(N*Nm ) 

-0.22 

-0.10 

1.77 

-0.03 

Draft  with 

:  (N*N )  or 

(N*Nm) 

0.07 

0.12 

-0.02 

A  z  with 

:  (N*Nm) 

-0.06 

-0.03 

Pitch  with 

:  (Nm*Nm) 

-0.10 

Corr.  Coefficient 

F  x  with 

% 

* 

-0.16 

-0.08 

0.65 

-0.06 

Draft  with 

* 

• 

0.16 

0.13 

-0.10 

Az  with 

m 

• 

-0.06 

-0.20 

Pitch  with 

• 

m 

-0.31 

SHEET  C 

VEL  = 

0.30 

(RUN  NO 

:  0301  8 

) 

Fx 

Draft 

Az 

Pitch 

Omega 

N 

m 

m/sec**2 

radi ans 

1 /sec**2 

MEAN 

11  .73 

C.QC20 

0.003 

0.0032 

-0.G02 

SIGMA 

2.76 

0.0003 

0.005 

0.0009 

0.044 

scaled  as 

forces  and 

moments 

N 

N 

N 

Nm 

Nm 

MEAN 

1 1  .73 

16.01 

0.73 

7.93 

-0.09 

SIGMA 

2.76 

2.45 

1.31 

2.12 

2.27 

Covariance 

of 

F  x  with 

:  (N*N  )  or  (N*Nm) 

1.19 

-1.30 

2.96 

-1.24 

Draft  with 

l  (N*N )  or  (N*Nm) 

-0.28 

1 .20 

-0.79 

Az  with 

:  (N*Nm) 

-0.83 

0.85 

Pitch  with 

:  (Nm*Nm) 

-3.57 

Corr.  Coefficient 

F  x  with 

* 

* 

0.18 

-0.36 

0.51 

-0.20 

D  raft  with 

♦ 

-0.09 

0.23 

-0.14 

Az  with 

* 

• 

-0.30 

0.29 

Pitch  with 

# 

• 

-0.74 

SHEET  D 

VEL  = 

0.03 

(RUN  NO 

:  04025 

) 

Fx 

Draft 

Az 

Fitch 

Omega 

N 

m 

m/sec**2 

radians 

1 /sec**2 

MEAN 

11.25 

0.0027 

-0.004 

0.0052 

0.003 

SI6MA 

3.70 

0.0002 

0.00? 

0.0015 

0.014 

scaled  as 

forces  and 

moments 

N 

N 

N 

Nm 

Nm 

MEAN 

11.25 

21.17 

-1 .00 

12.65 

0.16 

SIGMA 

3.70 

1  .2A 

0.47 

3.62 

0.74 

C  ovariance 

of 

Ex  with 

:  (N*N )  or  (N*Nm) 

1.08 

-0.80 

7.70 

-0.21 

Draft  with 

:  (N*N )  or  (N*Nm> 

-0.20 

2.68 

-0.07 

A  z  with 

r  (N*Nm) 

-0.95 

0.04 

Pitch  with 

:  (Nm*Nm) 

-0.6? 

Corr.  Coefficient 

F  x  with 

0.23 

-0.46 

0.58 

-0.08 

Draft  with 

-0.34 

0.60 

-0.08 

Az  with 

-0.56 

0.11 

Pitch  with 

-0.23 

SHEET  D 

VEL  = 

0.30 

(RUN  NO 

:  05025 

) 

Fx 

Draft 

Az 

Pitch 

Omega 

N 

m 

m/sec**? 

radians 

1 /se  c**  2 

MEAN 

15.16 

0.0024 

0.002 

0.0052 

-0.003 

SIGMA 

3.41 

C .0004 

0.009 

0.0013 

0 .060 

scaled  as 

forces  and 

moments 

N 

N 

N 

Nm 

Nm 

MEAN 

15.16 

19.11 

0.48 

12.80 

-0.13 

SIGMA 

3.41 

3.23 

2.21 

3.29 

3.10 

C  ovariance 

of 

F  x  with 

:  (N*N )  or  ( N*Nm ) 

4.14 

-1  .77 

5.59 

-2.01 

Draft  with 

:  (N*N)  or  (N+Nm) 

-2.27 

7.24 

-2.52 

A  z  with 

:  (N*Nm) 

-2.72 

3.82 

Pitch  with 

X  (Nm*Nm) 

-6.68 

Corr.  Coefficient 

F  x  with 

f 

» 

0.38 

-0.23 

0.50 

-0.19 

Draft  with 

• 

• 

-0.3? 

0.68 

-0.25 

A  z  with 

• 

* 

-0.37 

0.56 

Pitch  with 

• 

* 

-0.66 

SHEET  E 

VEL  =  0.03 

(RUN  NO 

:  0402 C 

) 

Fx 

Draft 

A  z 

Pitch 

Omega 

N 

m 

m/sec**2 

radians 

1  /sec**2 

MEAN 

10.46 

0.0015 

0.000 

0.0039 

0.000 

SIGMA 

3.34 

0.0001 

0.001 

0.00U5 

0.003 

scaled  as 

forces  and  moments 

N 

N 

N 

Nm 

Nm 

MEAN 

10.46 

12.09 

0.01 

9.60 

0.01 

SIGMA 

3.34 

0.68 

0.36 

1.26 

0.18 

C  ovar lance 

of 

Fx  with 

:  (N*N )  or  (N*Nm) 

-0.42 

-0.39 

1.65 

0.00 

Draft  with 

:  (N *N )  or  (N*Nm) 

0.02 

0.40 

-0.01 

A  2  with 

I  (N*Nm) 

0.02 

0.01 

Pitch  with  :  (hm*Nm)  -Q.03 


Com.  Coefficient 


Fx  with  : 

-0.18 

-0.32 

0.79 

0.00 

Draft  with  : 

0.07 

0.47 

-0.08 

A  z  wi th  : 

0.05 

0.11 

Pitch  with  : 

-0.15 

SHEET  E 

VEL  =  0.59 

(RUN  NO 

:  0502C 

) 

Fx 

Draf  t 

A  z 

Pitch 

Omega 

N 

m 

m/ sec**2 

radians 

1  /  sec**2 

MEAN 

18.63 

Q .001 1 

-0.001 

C . 00  56 

Q.001 

SIGMA 
scaled  as 

4.10 

forces  and  moments 

0.0002 

0.013 

0.0006 

0.056 

N 

N 

N 

Nm 

Nm 

MEAN 

18.63 

3.96 

-0.17 

13.78 

0.03 

SIGMA 

Covariance 

4.10 

of 

1.50 

3.18 

1.41 

2.88 

F  x  with 

I  (N*N)  or  (N*Nm) 

0.78 

0.36 

1.65 

0.36 

Draft  with  :  (h*N)  or  (N*Nm) 
Az  with  :  (N*Nm) 

Pitch  with  l  (Nm*Nm) 

Cor r .  Coef f i cient 

-2.16 

0.37 

-0.79 

-1  .24 
5.71 
-2.30 

Fx  with 

• 

• 

0.13 

0.03 

0.28 

0.03 

D  raft  with 
Az  with 
Pitch  with 

« 

« 

• 

• 

• 

-0.45 

0.18 

-0.18 

-0.25 

0.62 

-0.57 

SHEET  F 

VEL  = 

0.03 

(RUN  NO 

:  06028 

) 

F  x 

D  raf  t 

Az 

Fitch 

Omega 

N 

m 

m/sec**2 

radians 

1  /sec**2 

MEAN 

13.60 

G .  00  24 

-0.003 

C.0G58 

-0.002 

SIGMA 

3.84 

0.0Q01 

0.010 

0.0012 

0.009 

scaled  as 

forces  and 

moment  s 

N 

N 

N 

Nm 

Nm 

MEAN 

13. 6Q 

1  8  .Q1 

-0.75 

14.33 

-0.11 

SIGMA 

3.84 

0.90 

2.56 

2.98 

0.46 

Covariance 

of 

F  x  with 

:  (N*N )  or  <N*Nm) 

0.26 

-5.84 

1 .77 

0.06 

Draft  with 

I  ( N*N  )  or  (N*hm) 

0.29 

0.58 

-0.02 

Az  with 

I  (  N  *  N  m  ) 

-2.93 

0.01 

Pitch  with 

:  (Nm*Nm) 

-0.26 

Corr.  Coefficient 

F  x  with 

0.07 

-0.59 

0.16 

0.04 

Draft  with 

0.11 

0.20 

-0.04 

Az  with 

-0.39 

0.01 

Pitch  with 

-0.19 

SHEET  F 

VEL  = 

0.60 

(RUN  NO 

:  07029 

) 

Fx 

Draft 

Az 

Pitch 

Omega 

N 

m 

m/sec**2 

radians 

1  / sec**2 

MEAN 

26.54 

0.0025 

-0.010 

0.0085 

-0.002 

SIGMA 

7.25 
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APPENDIX  B 


Time  Series  of  the  Measured  Quantities 
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APPENDIX  C 

Power  Spectra,  Cross-Spectra  and  Coherance-Functions 


Figure  Cl:  Power  Spectra 

a)  Examples  of  low  velocity  runs,  VQ  =  .03  m/sec 

b)  Runs  at  medium  velocity,  V0  =  0.22  m/sec  and  VQ  =  0.3  m/sec 

c)  Runs  at  high  velocity,  VQ  =  0.6  m/sec 

d)  Runs  with  Sheet  G,  VQ  =  0.03  m/sec,  VQ  =  0.22  m/sec 

Figure  C2:  Cross  Spectra  and  Coherance  functions  of  z(t) 

with  z(t),  0 ( t )  and  0 ( t ) 

a)  Examples  of  low  velocity  runs,  VQ  =  0.03  m/sec 

b)  Runs  at  medium  velocity,  VQ  =  .22  m/sec  and  V0  =  0.3  m/sec 

c)  Runs  at  high  velocity,  VQ  =  0.6  m/sec 

d)  Runs  with  sheet  G,  V0  =  0.03  m/sec  V0  =  0.22  m/sec 


Figure  C3:  Cross  Spectra  and  Coherence  Functions  of  Fx(t)  with  z(t), 

Fx(t)  with  z(t) 

a)  Examples  of  low  velocity  runs,  VQ  =  0.03  m/sec 

b)  Runs  at  medium  velocity,  V0  =  0.22  m/sec  and  VQ  =  0.3  m/sec 

c)  Runs  at  high  velocity,  V0  =  0.6  m/sec 

d)  Runs  with  sheet  G,  VQ  =  0.03  m/sec  and  VQ  =  0.22  m/sec 

Figure  C4:  Cross-Spectra  and  Coherence  Functions  of  Fx(t) 

with  © ( t)  and  F  (t)  with  6(t) 

a)  Examples  of  lew  velocity  runs,  VQ  =  0.03  m/sec 

b)  Runs  at  medium  velocity,  VQ  =  0.22  m/sec  and  VQ  =  0.3  nv/sec 

c)  Runs  at  high  velocity,  V0  =  0.6  m/sec 

d)  Runs  with  sheet  G,  VQ  =  0.03  m/sec  and  VQ  =  0.22  nv/sec 
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Figure  C2:  Cross-Spectra  and  Coherance  Functions  of  Draft  z(t)  with 

Vertical  Acceleration  z(t),  and  Pitch  0(t)  with  Angular 
Acceleration  6 ( t ) 
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Figure  C3:  Cross-Spectra  and  Coherance  Functions  of  Fx(t)  with 

draft  z(t) ,  and  Fx(t)  with  Vertical  Acceleration  z(t) 
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Figure  C4:  Cross-Spectra  and  Coherance  Functions  of  Fv(t)  with 

Pitch  0(t)  and  Fx(t)  with  Angular  Acceleration,  e(t) 
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